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1)  FOREWORD 


This  final  report  covers  the  work  performed  by  Vincent  G.  Dominic  over  the  period  of 
August  1989  to  August  1992  under  the  auspices  of  the  Joint  Services  Electronics 
Program  (JSEP)  graduate  fellowship.  During  this  fellowship,  we  investigated  two 
distinct  problems.  First,  we  studied  short-pulse  interactions  in  photorefractive  crystals 
and  the  use  of  a  photorefractive  crystal  to  measure  the  coherence  length  of  laser  pulses. 
This  work  has  been  published  in  two  separate  articles b2  and  reprints  are  included  in 
Appendix  B  of  this  report.  Second,  we  studied  second-harmonic  generation  in  optical 
fibers  and  glass  optical  fiber  preforms.  This  topic  will  be  the  main  emphasis  of  this  final 
report. 

2)  INTRODUCTION 

During  my  JSEP  fellowship  tenure  we  studied  two  completely  separate  problems.  The 
first  problem  was  how  short  laser  pulses  interact  inside  a  photorefractive  crystal,  such 
as  barium  titanate  (BaTiOa),  and  in  particular  the  use  of  photorefractive  crystals  to 
shape  the  temporal  profile  of  short  pulses  of  light.  We  also  developed  a  nice  technique 
to  determine  the  coherence  length  of  mode-locked  laser  pulses  in  real  time.  Our 
coherence-length  measurement  technique  works  well  for  pulses  having  coherence  times 
ranging  from  100  ps  down  to  about  1  ps,  as  discussed  in  Ref.  1.  We  also  proposed  a 
theoretical  model  to  describe  short  pulse  interactions  inside  photorefractive  crystals.^ 
Although  I  have  not  been  directly  involved  in  the  continuation  of  this  work,  recent 
experiments  have  born  out  the  validity  of  our  theory.^ 

The  second  and  primary  research  work  sponsored  by  this  JSEP  fellowship  is  to 
understand  second-harmonic  generation  in  glass.  This  effect  was  first  observed  in 
optical  fibers'^  but  also  appears  in  bulk  glass  samples.5-7  This  effect  is  interesting 
because,  at  first  glance,  it  appears  to  be  impossible.  Glass  is  centrosymmetric,  and  thus 
second-order  nonlinear  effects  are  forbidden  in  the  dipole  approximation.®  However,  in 
1986  Osterberg  and  Margulis'*  observed  that  if  they  launched  intense  1.064  pm  laser 
pulses  from  a  Q-switched  Nd:YAG  laser  into  an  optical  fiber,  there  was  a  slow  buildup 
of  green  532  nm  light  produced  by  the  fiber.  The  buildup  time  was  several  hours. 
Stolen  and  Tom^  later  showed  that  if  some  externally  generated  532  nm  light  is 
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launched  into  the  fiber  along  with  the  fundamental  light,  then  the  buildup  time  is 
reduced  from  hours  to  minutes. 

What  happens  inside  the  glass  to  remove  the  centrosymmetry  of  the  material?  It  is  now 
agreed  that  a  DC  electric  field  is  induced  inside  the  glass. This  field  produces  a 
second-order  nonlinearity  via  the  third-order  susceptibility  by  electric  field  induced 
second-harmonic  generation  (EFISH).  Thus  the  induced  nonlinearity  can  be  written: 

X^-\-2(o-,co.(o)  -  3  x^^\-2co\(o,(o,0)  E^  (l) 

where  is  the  known  third-order  susceptibility  tensor  of  glass,  is  the  induced  DC 
electric  field,  and  is  the  resulting  second-order  nonlinearity.  Here  we  follow  the 
notation  of  Butcher  and  Cotter. ’3  The  polarization  properties  of  the  induced 
nonlinearity  have  been  shown  to  follow  the  prediction  of  Eq.  (1)  in  optical  fibers^i^  and 
in  bulk  glass  samples,^ 3  although  the  spatial  patterns  of  light  generated  according  to 
Eq.  (1)  can  be  more  easily  seen  in  bulk  glasses  than  in  fibers,  and  this  has  been  one  of  the 
main  thrusts  of  our  research.  We  have  found'^  that  Eq.  (1)  is  fully  born  out,  and  that  it 
is  possible  to  use  its  polarization  properties  to  map  out  the  spatial  profile  of  the  induced 
electric  field.  In  all  of  our  work  we  have  found  it  safe  to  assume  that  x^^^  in  Eq.  (1) 
obeys  Kleinman  symmetry,  as  was  also  found  to  be  true  in  Ref.  10. 

What  happens  inside  the  glass  to  cause  the  formation  of  this  DC  electric  field? 
Although  this  question  has  not  been  fully  answered,  and  is  the  subject  of  ongoing 
debate,  a  majority  of  researchers  now  believe  that  the  interference  of  two  or  more  multi¬ 
photon  ionization  channels  gives  rise  to  a  spatially  anisotropic  photoelectron 
current.^'^'l^  This  photogalvanic  current  continues  until  the  DC  electric  field  it  produces 
becomes  large  enough  to  inhibit  any  further  separation  of  electrical  charges.  Figure  1 
shows  two  different  models  which  lead  to  this  photogalvanic  current.  Model  A  shows 
the  interference  of  one-photon  and  two-photon  excitation  of  an  electron,  and  the 
anisotropy  in  the  direction  of  the  ejected  electron  has  been  convincingly  demonstrated 
in  gases. Model  B  is  believed  to  be  more  likely  in  our  glass  samples,  in  which  there 
is  significant  absorption  at  the  energy  equal  to  four  infrared  (1.064  pm)  photons.  This 
absorption  band  in  glass  is  associated  with  its  germanium  dopant,  and  germanium 
doping  is  necessary  to  observe  second-harmonic  generation  in  fused  silica.  Because 
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Model  A  Model  B 


Figure  1)  Energy  level  model  associated  with  anisotropic  photoelectron  emission  from  germanium- 
doped  fused  silica.  Model  A  is  the  simple  picture,  however  there  is  very  little  absorption  at 
X=532  nrn  in  the  glass  materials  generally  used  for  these  experiments.  In  model  B  there  are 
three  separate  ionization  channels  that  interfere.  The  defect  level  is  associated  with  the 
known  >.=240  nm  absorption  band  in  germanium-doped  glass. 

model  B  applies  to  germanium-doped  glass  systems  (which  includes  most 
communication  optical  fibers  and  preforms),  we  expect  the  second-harmonic  generation 
effect  to  be  strongly  intensity  dependent,  due  to  the  multi-photon  nature  of  the 
ionization.  Consequently,  the  second  major  thrust  of  the  research  performed  under  this 
JSEP  fellowship  was  to  measure  the  intensity  dependence  of  this  effect  to  see  if  it  fit 
model  B.  We  found  that  the  intensity  dependence  was  even  stronger  than  expected 
from  model  B,i^  and  that  the  photo-ejection  of  electrons  is  more  complicated  than 
predicted  by  this  simple  model.  In  particular  we  found  evidence  that  other  energy 
levels  (metastable)  of  the  glass  system  are  involved.  Our  paper  describing  this  work^-^ 
has  been  accepted  for  publication  in  Optics  Letters  (see  Appendix  A). 

3)  SUMMARY  OF  RESULTS 

There  are  three  important  conclusions  from  our  work: 

*  Experiments  involving  second-harmonic  generation  in  glass  are  simpler  and  easier 
to  interpret  if  they  are  performed  in  bulk  glass  samples  rather  than  in  optical  fibers. 
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We  used  optical  fiber  preform  samples  made  from  the  same  glass  as  the 
communication  fibers  used  by  most  researchers. 

*  The  creation  of  the  second-order  nonlinearity  inside  the  glass  can  be  completely 
explained  by  invoking  a  DC  electric  field  Thus,  the  following  ecjuation  explains 
the  induced  nonlinearity: 

X^^\-2co-,co,co)  =  ?>  x^'‘^i-2(o;co,a},0)  E^. 

*  The  mechanism  that  induces  the  DC  electric  field  definitely  involves  multi-photon 
interactions.  This  is  demonstrated  by  our  measurement  of  the  intensity  dependence 
of  the  rate-of-formation  of  the  DC  field.  It  seems  likely  that  interference  of  multi¬ 
photon  ionization  channels  in  the  presence  of  a  laser  beam  (Eq))  and  its  second 
harmonic  (E2(o)  produces  an  anisotropic  distribution  of  ejected  electrons,  which  can 
be  thought  of  as  a  photogalvanic  current.  This  current  builds  up  the  DC  electric  field 

4- 

4)  ADVANTAGES  OF  BULK  GLASS  SAMPLES  OVER  OPTICAL  FIBERS 

We  began  our  studies  of  second-harmonic  generation  using  optical  fibers.  However, 
after  struggling  for  about  a  year  and  a  half  to  obtain  reproducible  results,  we  switched 
to  experiments  in  bulk  glass  samples.  Our  bulk  glass  sample  is  an  optical  fiber  preform. 
Optical  fiber  preforms  are  carefully  grown  cylindrical  glass  rods  which  are  later  heated 
to  the  melting  point  and  drawn  into  optical  fibers.  For  the  purposes  of  our  experiments, 
the  bulk  glass  preform  is  equivalent  to  an  optical  fiber.  There  are  minor  differences  in 
the  glass  structure  of  a  fiber  versus  a  preform,  such  as  the  number  density  of  certain 
kinds  of  defects,^^  however  we  don't  expect  this  to  adversely  affect  the  applicability  of 
our  results.  There  are  definite  advantages  to  using  bulk  glass  samples  over  fibers  for 
these  experiments.  These  advantages  made  it  possible  to  carry  out  reproducible,  well- 
controlled  experiments.  I  list  below  the  major  advantages. 

Elimination  of  modal  overlap  problems 

Only  certain  spatial  distributions  of  light  (modes)  can  propagate  in  an  optical  fiber,  and 
this  causes  problems.  For  example,  in  our  experiments  we  need  to  inject  two  different 
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color  light  beams,  one  at  frequency  co  and  one  at  frequency  2a),  and  the  allowed  modal 
shapes  for  the  two  different  colors  will  be  different  in  a  fiber.  Even  if  a  fiber  is  single 
mode  at  to,  in  general  it  will  be  multi-mode  at  2co.  This  means  that  the  spatial  overlap  of 
the  two  optical  fields  (E^j  and  E2co)  ^^’bl  change  if  the  coupling  into  the  optical  fiber 
changes  or  if  the  fiber  is  perturbed  by  the  slightest  vibration  or  bend.  Consequently,  the 
spatial  overlap  depends  on  the  exact  launching  conditions  into  the  fiber,  and  can  change 
from  one  experiment  to  the  next,  or  even  during  a  single  experiment.  Because  the 
modal  overlap  of  the  launched  E(j  and  E2a)  fields  is  difficult  to  control,  it  makes  it  quite 
difficult  to  obtain  consistent,  reproducible  experimental  results.  There  is  one  remedy 
for  this  situation:  use  a  fiber  that  is  single  mode  at  both  O)  and  2co.  However,  in  this 
case  the  fundamental  (E^o)  wavelength  is  only  weakly  guided  and  the  majority  of  its 
energy  actually  propagates  in  the  cladding  region  of  the  fiber.  This  loss  of  power  in  the 
core  of  the  fiber  greatly  reduces  the  rate  at  which  electrons  are  ejected  in  the  fiber,  and  it 
also  makes  it  tricky  to  determine  the  amount  of  power  actually  contained  in  the  core, 
because  power  measurements  are  easily  polluted  by  light  propagating  down  the  fiber  in 
so-called  "cladding  modes."  In  the  case  of  such  weak  guiding,  the  fundamental  optical 
field  is  also  quite  susceptible  to  bending-induced  out-coupling,  so  that  small  vibrations 
cause  the  light  to  leave  the  optical  fiber  core  and  radiate  away.  Because  of  the 
complications  caused  by  the  use  of  a  fiber  that  is  single  mode  at  both  co  and  2(o,  only  one 
experiment  has  been  reported  to  date  in  such  a  fiber. B) 

Another  restriction  imposed  by  the  discrete  modes  of  an  optical  fiber,  and  which  is 
absent  in  a  bulk  sample,  is  the  spatial  shape  of  the  generated  second-harmonic  signal.  I 
will  report  later  in  this  document  how  the  spatial  shape  of  the  generated  second- 
harmonic  beam  can  give  very  useful  information  about  the  shape  of  the  internal  DC 
electric  field.  Such  experiments  would  be  impossible  in  an  optical  fiber,  for  not  only  is 
the  spatial  shape  of  the  generated  signal  restricted  to  the  allowed  modes  of  the  fiber,  but 
also  the  fiber  geometry  does  not  allow  one  to  probe  the  induced  field  at  any  other 
point  except  at  the  center  of  the  fiber.  This  means  that  in  a  fiber  one  cannot  map  the 
shape  of  the  induced  DC  field  by  translating  a  probe  beam,  as  we  can  in  bulk  glass 
samples,  as  we  describe  in  Section  5  below. 
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Elimination  of  de-phasing  problems 

Since  we  must  focus  rather  tightly  into  a  bulk  glass  sample  to  achieve  the  same  optical 
intensities  as  exist  in  optical  fibers,  we  have  correspondingly  short  interaction  lengths. 
This  means  that  the  measured  second-harmonic  signal  will  be  weak.  Howev'er,  the 
reduced  interaction  length  (200  pm  versus  20  cm)  means  that  the  interaction  will  be 
about  100  times  less  sensitive  to  dephasing  effects.  There  are  several  important 
dephasing  mechanisms,  including  self-phase  modulation,  cross-phase  modulation,  and 
temperature-induced  changes  in  the  refractive  index;  all  of  these  effects  were  found  to 
be  present  in  our  experiments  with  optical  fibers  and  absent  in  our  experiments  with 
bulk  glass. 

Various  other  advantages 

Another  major  advantage  of  bulk  samples  is  that,  once  they  are  cut  and  polished,  it  is 
possible  to  perform  many  independent  experiments  in  different  regions  of  the  sample 
before  it  is  necessary  to  switch  to  a  new  sample.  This  is  not  true  of  optical  fibers,  which 
must  be  discarded  after  each  experiment.  This  means  that  exorbitant  amounts  of  time 
are  spent  on  fiber  preparation  and  alignment.  Also,  because  the  required  intensities  are 
very  near  the  damage  threshold  of  fused  silica,  the  input  cleave  of  a  fiber  is  frequently 
damaged.  In  a  bulk  sample,  if  a  region  becomes  damaged  one  simply  translates  the 
sample  away  from  the  damage  spot.  Also,  at  low  intensities  the  rate-of-formation 
becomes  very  slow  (hours),  and  it  is  difficult  to  maintain  coupling  stability  into  a  fiber 
for  such  extended  times. 

In  all  fairness,  there  are  two  drawbacks  to  using  bulk  samples  which  should  be  pointed 
out.  First,  their  relatively  short  interaction  length  limits  the  generated  second-harmonic 
signal  to  quite  small  values,  such'  that  sensitive  detectors  (PMTs)  and  electronics  must 
be  used.  This  is  not  a  major  difficulty.  Second,  precise  spatial  overlap  of  the 
fundamental  and  second-harmonic  seeding  beams  requires  that  special  optics  be  used 
to  focus  the  laser  into  the  sample.  This  is  a  rather  tricky  experimental  point  which  is  not 
widely  appreciated. 
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5)  ELECTRIC-FIELD  INDUCED  SECOND-HARMONIC  GENERATION 


Mizrahi  ct  al.  showed  in  1990  that  the  polarization  properties  of  the  induced 
nonlinearity  could  be  completely  explained  by  invoking  the  known  nonlinear 
susceptibility  tensor  of  glass  and  a  single  cartesian  component  of  the  induced  DC 
electrical  fieldd^  Their  experiments  were  carried  out  in  a  fiber  that  was  single  mode  at 
both  the  fundamental  (co)  and  the  second-harmonic  frequency  (2a)).  They  seeded  the 
sample  with  both  beams  polarized  vertically,  which  should  cause  charges  to  be 
separated  vertically  by  the  photogalvanic  current.  To  explain  their  observations, 
Mizrahi  et  a},  needed  only  to  invoke  the  vertical  component  of  the  DC  electric  field. 
There  are  two  reasons  for  this;  1)  They  can  only  launch  light  into  the  core  of  the  optical 
fiber,  and  since  this  is  the  point  of  symmetry  for  the  vertical  charge  separation,  the  field 
is  mostly  vertically  oriented  in  this  region.  2)  The  symmetry  of  the  horizontal 
component  of  the  DC  electric  field  is  such  that  it  will  generate  a  four  lobed  mode  that  is 
not  an  allowed  propagation  mode  in  this  fiber.  (Thus  the  single-mode  nature  of  the 
fiber  prevents  determination  of  the  shape  of  the  DC  electric  field.) 

We  have  extended  the  work  of  Mizrahi  t’t  al.  by  performing  similar  experiments  in  bulk 
glass  samples.'-  In  the  bulk  sample  we  are  free  to  probe  the  induced  field  at  whatever 
point  we  choose  merely  by  translating  a  probe  beam.  Also,  the  shape  of  the  generated 
second-harmonic  signal  is  not  constrained.  We  performed  two  independent 
experiments  to  determine  the  shape  of  the  induced  electrical  field.  In  the  first 
experiment  we  seeded  the  bulk  sample  with  an  actively  mode-locked  (76  MHz,  100  ps 
FWHM),  Q-switched  (1  kHz,  300  ns  FWHM)  laser  beam  which  was  doubled  in  a 
temperature-tuned,  non-critically  phase-matched  crystal  of  LBO  (LiB03).  The  average 
power  in  the  infrared  (1.064  pm)  was  typically  4  W  and  in  the  green  (0.532  pm)  was 
~2  mW.  Both  colors  were  focused  to  the  same  location  inside  the  sample  with  a  5X  lens 
designed  specifically  to  be  achromatic  at  these  two  wavelengths.  The  infrared  and 
green  beams  had  1/e-  of  intensity  radii  of  ~9  pm  and  ~7.5  pm,  respectively.  The  sample 
was  seeded  for  3  minutes,  during  which  time  the  generated  second-harmonic  signal  was 
monitored  until  it  reached  a  nearly  saturated  level.  Both  laser  beams  were  vertically 
polarized  during  seeding.  After  3  minutes,  we  blocked  the  seeding  green  beam  and 
decreased  the  infrared  power  so  that  no  further  change  was  induced  inside  the  glass. 
We  inserted  a  X/l  waveplate  (@  1.064  pm)  immediately  before  the  focusing  achromat 
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Figure  2)  Experimental  setup  to  map  out  the  transverse  dependence  of  the  induced  DC  electric  field 
inside  a  glass  optical  fiber  preform. 

and  an  analyzer  after  the  collecting  objective  (see  Fig.  2),  so  that  we  could  control  the 
input  infrared  reading  polarization  and  the  detected  signal  polarization.  The  generated 
signal  was  proportional  to: 

^2(0  I  ^Inalyzer  ^  co  ^  CO  ^dc)\  (2) 

where  is  the  infrared  probing  beam,  is  the  DC  electric  field  and  e^n^iyzer  is  the  unit 
vector  describing  the  orientation  of  the  analyzer.  Because  of  the  symmetry  of  iri 
glass,  if  we  vertically  polarize  the  probe  beam  and  we  set  the  analyzer  to  vertical  then 
only  the  vertical  component  of  can  contribute  to  the  measured  signal.  Thus,  if  we 

I'  ~  1^ 

move  the  probe  beam  around  we  can  map  the  strength  of  as  a  function  of 

transverse  position.  We  can  also  perform  the  complementary  experiment  in  which  the 
probe  beam  is  horizontally  polarized  and  the  analyzer  is  horizontally  oriented.  This 
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maps  the  horizontal  component  |i  E,^  j  ot'  tlie  DC  electric  field.  We  take  the  square  root 

of  such  data  to  obtain  a  map  of  the  two  transverse  components  of  the  DC  field.  To  turn 
this  information  into  a  \ector  map  of  the  field  we  need  to  perform  a  separate 
experiment  to  determine  the  relative  sign  between  the  \'ertical  and  horizontal  field 
components.  We  determine  the  relative  phase  between  the  transverse  field  components 
by  rotating  the  input  probe  polarization  with  a  fixed  analyzer  (either  vertical  or 
horizontal).  The  shape  of  the  signal  as  the  input  polarization  is  rotated  determines 
whether  the  product  E"l E'.'"'  is  positive  or  negative 


There  is  a  complication  to  this  experiment.  Since  the  probing  beam  was  not  infinitely 

I  —  j  1  -•  i'' 

and  over  the  area  of  the 

probing  beam,  which  smears  out  the  information.  Also,  since  the  sign  of  the  term 
E^'‘  E'//‘  may  be  different  in  different  regions  of  the  probed  area,  we  cannot  uniquely 
assign  either  a  positic’e  or  a  negative  sign  to  this  term.  A  simple  way  to  overcome  this 
difficulty  is  to  allow  the  £^."‘  EV.’’’  term  to  be  complex  (neither  fully  positive  or 
negative).  This  works  wonderfully.  In  this  case,  the  cross-product  term  may 

be  assigned  a  relative  phase  o  which  enters  the  prediction  for  the  measured  signal  only 
as  cos(o].  If  0>:t/2  then  the  cross-product  term  is  more  negative  than  positive,  and  the 
opposite  is  true  if  o<:t/2.  In  this  way  we  can  use  our  experiments  with  the  rotating 
probe  polarization  to  determine  the  relative  sign  between  the  vertical  and  horizontal  DC 
field  components  in  the  glass  caused  by  its  previous  irradiation  by  light  beams. 


Here  are  our  results.  We  first  assume  a  shape  for  the  transverse  distribution  of  charge. 
This  charge  distribution  is  assumed  to  be: 


p{x,\) 


-Po  vexp 


.r  +  1 


(5  pw)- 


(3) 


which  is  simply  the  y-derivative  of  a  cylindrically  symmetric  gaussian.  This  leads  to  the 
field  distribution  shown  in  Fig.  3.  We  wrote  a  simple  computer  program  to  predict  how 
the  measured  second-harmonic  signal  should  vary  as  the  probe  beam  location  is  varied. 
This  program  assumed  that  the  DC  electric  field  was  as  shown  in  Fig.  3,  and  it  took  into 
account  the  finite  size  of  the  probing  beam.  Figures  4  and  5  show  a  comparison  of  the 
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Figure  3)  This  is  the  DC  electric  field  of  the  charge  distribution  given  by  Eq.  (3)  in  the  text.  This  electric 
field  is  used  to  predict  the  results  of  two  separate  experiments  in  our  bulk  glass  samples. 
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Predicted 


Observed 


Horizontal  Position  ( \xm  ) 


Horizontal  Position  ( jim  ) 


Figure  4)  Comparison  between  the  theoretically  predicted  and  the  experimentally  observed  second- 
harmonic  signal  versus  position  of  the  probing  beam.  The  polarization  of  the  probing  beam 
is  vertical  and  the  analyzer  passe's  only  vertically  polarized  light.  These  graphs  show  contour 
plots  of  equal  signal  strength.  The  strength  of  the  contour  values  increase  monotonically 
towards  the  center  of  the  graphs.  The  wings  on  the  top  and  bottom  of  the  contour  maps 
correspond  to  regions  where  the  vertical  component  of  the  DC  electric  field  switches  sign. 
The  theoretical  map  was  generated  using  the  DC  electric  field  shown  in  Fig.  3  and  the  known 
9  pm  radius  of  the  probing  beam. 
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Vertical  Position  ( jim  ) 


Predicted  Observed 


Horizontal  Position  ( |j.m  )  Horizontal  Position  ( |j.m  ) 


Figure  5)  Same  as  Fig.  4  except  now  the  probing  light  is  polarized  horizontally  and  the  analyzer  is  set 
to  pass  horizontally  polarized  light.  In  this  case  the  contour  plots  show  a  four-lobed  pattern. 
(With  a  much  smaller  probing  beam  the  signal  would  drop  to  zero  between  each  of  the  four 
peaks.)  The  slight  asymmetry  in  the  experimental  data  is  attributed  to  imperfect  overlap  of 
the  fundamental  and  second-harmonic  seeding  beams.  The  theoretical  map  was  generated 
using  the  DC  electric  field  shown  in  Fig.  3  and  the  known  9  pm  radius  of  the  probing  beam. 
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predicted  and  observed  dependence  of  the  strength  of  the  generated  second-harmonic 
signal  as  the  position  of  the  probing  beam  is  translated.  Figure  4  shows  the  case  where 
the  probe  polarization  is  vertical,  the  analyzer  orientation  is  vertical,  and  the  preform 
sample  was  translated  side  to  side  and  up  &  down  to  map  out  the  spatial  dependence. 
Figure  5  gives  the  corresponding  comparison  for  the  horizontal  input 
polarization/ horizontal  analyzer  case.  The  comparison  was  quite  good.  However  there 
was  a  significant  amount  of  smearing  of  the  spatial  information  because  of  the  large 
probe  beam  size.  This  is  especially  obvious  in  Fig.  5.  When  we  rotated  the  probe 
polarization  with  a  fixed  analyzer  to  determine  the  sign  of  Ej'.'''  in  each  of  the  four 
quadrants,  we  confirmed  that  the  DC  field  components  have  the  relative  sign  shown  in 
Fig.  3.  In  brief,  then,  these  experiments  confirm  that  the  induced  DC  electric  field  looks 
like  the  field  shown  in  Fig.  3.  The  asvmmetrv  in  the  observed  data  of  Fig.  5  is  thought 
to  arise  from  the  slightly  imperfect  spatial  overlap  of  the  infrared  and  green  seeding 
beams. 

We  also  performed  an  independent  experiment  to  confirm  that  the  DC  field  had  the 
assumed  shape.  We  observed  the  spatial  shape  of  the  generated  second-harmonic 
signal  as  the  probe  polarization  is  rotated,  as  first  suggested  by  Driscoll  and  Lawandy.^^ 
We  observed  the  signal's  mode  shape  using  either  a  vertical  or  horizontal  analyzer,  as 
shown  in  Fig.  6.  For  this  experiment  we  would  prefer  to  have  an  infinitely  big  probing 
beam.  However,  since  this  is  not  feasible  we  can  use  the  same  computer  program 
discussed  previously  to  take  into  account  the  finite  size  of  the  reading  beam.  The 
program  predicted  how  the  far-field  mode  shape  should  vary  with  the  reading 
polarization. 

We  made  a  video  which  shows  a  side-by-side  comparison  between  our  computer 
prediction  of  the  far-field  mode  shape  and  the  actually  observed  mode  shape  dance. 
We  animated  the  computer  prediction  and  recorded  the  computer  screen.  We  then 
used  the  setup  shown  in  Fig.  6  to  film  the  observed  signal  on  the  CCD  camera.  In  both 
cases,  using  either  a  vertically  or  a  horizontally  oriented  analyzer,  the  agreement 
between  the  prediction  and  the  observation  was  quite  good.  This  video  premiered  at 
the  Optical  Society  of  America  Annual  Meeting  in  September  1992,  and  is  available 
upon  request. 
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Figure  6)  Experimental  setup  to  observe  the  far-field  spatial  mode  pattern  of  the  generated  second- 
harmonic  signal.  The  analyzer  is  set  to  pass  either  vertically  or  horizontally  polarized  light, 
and  the  input  polarization  is  continuously  rotated  using  a  X/2  plate  mounted  on  a  motorized 
rotation  stage.  The  mode  pattern  is  viewed  with  a  CCD  camera  and  a  television  monitor. 

Two  independent  experiments  confirm  that  the  induced  electric  field  looks  as  shown  in 
Fig.  3.  We  consider  this  a  complete  verification  of  the  predictions  of  Eq.  (1)  above.  We 
now  have  no  doubt  that  second-harmonic  generation  in  glasses  is  caused  by  a  DC 
electric  field  inside  the  material. 

6)  MULTI-PHOTON  NATURE  OF  DC  FIELD  FORMATION  MECHANISM 

We  have  performed  extensive  measurements  of  the  rate-of-formation  of  the  DC  electric 
field  in  a  germanium-doped  optical  fiber  preform.  We  presented  this  work  at  CLEO 
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1992  (Ref.  6)  and  vve  expect  it  to  be  published  in  Optics  Letters  soon  (Ref.  7).  Since  our 
submission  to  Optics  Letters  contains  all  the  pertinent  details,  I  include  a  pre-print  of 
this  paper  here  as  an  appendix. 
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Growth  rate  of  second-harmonic  generation  in  glass 


V'ince  Dominic  and  Jack  Feinbcrg 
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Wc  study  iho  growth  rate  of  tho  light-induced  second-order  nonlinearity  in  germanium-doped  optical  fiber 
preforms.  We  seed  the  glass  with  bodi  infrared  light  and  green  light  (Nd:YAG  and  doubled  Nd:YAG)  to 
create  the  second-order  nonlinearity,  and  measure  its  formation  rate  while  varying  the  intensity  of  either 
the  fundamental  or  second-harmonic  seeding  beams  We  find  that  the  formation  rale  varies  as  a  power  law 
of  the  intensities,  but  with  an  exponent  larger  than  predicted  by  recent  models. 


Because  ihe  core  of  an  optical  fiber  is  inacroscopically 
cenlrosymmetric,  iwo-wave  sum  frequency  generation  fanil 
in  particular  second-harmonic  generation)  is  not  expected  to 
occur.  However,  in  1987  Osterberg  and  Margulis  reported 
rather  efficient  (5%)  second-harmonic  generation  of  1.064 
pm  light  inside  a  germanium-doped  optical  fiber  that  had 
been  illuminated  by  intense  infrared  light  for  many  hours.' 
Stolen  and  Tom^  showed  that  by  injecting  both  infrared  and 
green  light  into  the  fiber,  the  preparation  time  of  the  fiber 
could  be  shortened  from  many  hours  to  just  a  few  minutes. 

What  breaks  the  inversion  symmetry  of  the  glass? 
Recent  experiments  indicate  that  the  incident  light 
gradually  produces  a  spatially  periodic^  DC  electric  field  in 
the  glass.**-^  This  DC  electric  field  induces  a  second-order 
nonlinearity: 

x'^^\-2(o\(o,(o)  =  fl) 

Phase-matching  can  be  achieved  if  the  DC  field  (E<,c) 
periodically  changes  direction  along  the  fiber  length.  Many 
theories  have  been  postulated  to  explain  the  origin  of  this 
spatially  periodic  DC  electric  ficld.’'^'  ’  However,  these 
different  models  predict  different  dependencies  of  the  growth 
rate  of  the  DC  field  on  the  intensity  of  the  injected  optical 
beams.  To  discriminate  among  these  models  we  measured 
the  dependence  of  the  growth  rate  of  the  .second-harmonic 
signal  on  the  intensity  of  the  injected  light  beams. 

We  found  it  simpler  to  perform  our  experiments  in  glass 
fiber  preforms  instead  of  in  the  fibers  themselves.  Our 
preform  sample  had  an  outer  diameter  of  12.5  mm  and  a 
3.4  mm  diameter  germanium-doped  core.'®  The 
germanium  concentration  decreased  abruptly  at  the  core¬ 
cladding  interface  and  then  fell  exponentially  with  distance, 
with  occasional  rings  of  higher  index.  These  rings  result 
from  the  fabrication  process,  and  they  aberrate  optical 
beams  transmitted  through  them.  We  performed  our 
experiments  far  enough  away  from  the  core  to  avoid  these 
aberrations.  We  worked  in  the  cladding  region  .so  that  we 
could  fit  many  more  experiments  into  one  rotation  of  the 
sample.  We  verified  that  the  dependence  of  the  growth  rate 
on  intensity  was  the  same  in  the  core  region  as  m  the 
cladding  region  of  the  preform.  Before  starting  each 
experiment  we  rouited  the  sample  to  a  new,  unused  portion 
of  the  preform.  Becau.se  the  doping  was  cylindrically 
symmetric,  the  germanium  density  remained  constant  from 
experiment  to  experiment. 


We  used  a  seeding  technique’  in  which  we  focused  both 
the  fundamental  and  some  externally  generated  second- 
harmonic  light  into  the  sample.  (In  the  absence  of  a 
second-harmonic  seeding  beam  no  second-harmonic 
generation  was  observed  from  the  preform.)  A  mode-locked 
and  Q-switched  Nd:YAG  laser  [Coherent  Antares]  had  a 
Q-swiich  repetition  rate  of  1  kHz,  with  -23  mode-locked 
pulses  under  the  300-ns  Q-switch  pulse  envelope.  Its 
1.064  pm  light  was  doubled  in  a  non-crilically  phase- 
matched  LBO  crystal  to  produce  second-harmonic  light 
colinear  with  the  fundamental  light.  We  controlled  the 
power  levels  of  the  fundamental  and  second-harmonic 
seeding  beams  independently  with  two  harmonic 
waveplates.  The  paths  of  our  seeding  fundamenul  (Iq,)  and 
sccond-hturnonic  (hJ  beams  were  never  separated  to  ensure 
that  their  relative  phase  remained  stable  during  seeding." 
Wc  checked  this  phase  stability  and  found  that  it  exceeded 
four  hours. 

A  .spinning  wheel  (-5  Hz)  in  front  of  the  sample 
periodically  and  momentarily  blocked  the  green  seeding 
beam  (with  both  a  cold  mirror  and  an  RG780  absorption 
filter),  so  that  any  green  light  generated  in  the  sample  could 
be  measured.  Light  rcllecied  from  the  cold  mirror  was  used 
to  trigger  the  data  acquisition  system.  The  power  in  each 
incident  infrared  and  green  pulse  was  monitored  with 
photodiodes,  and  the  generated  green  light  was  monitored 
with  a  photomultiplier.  All  of  these  detectors  were 
carefully  checked  for  linearity  over  theu"  range  of  detected 
intensities.  A  25  mm  focal  length  achromatic  lens  focused 
the  infrared  and  the  green  to  the  same  location  within  the 
glass  preform.  The  infrared  and  green  beams  were  measured 
to  be  9.7  pm  and  7.4  pm  (radius  at  1/e^  of  intensity), 
respectively,  at  the  focus. 

We  acquired  two  types  of  data  sets:  1)  we  held  the 
infrared  power  constant  and  varied  the  green  power  level,  or 
2)  we  held  the  green  constant  and  varied  the  infrared  power 
level.  For  both  types  of  experiments,  we  measured  the 
early  time  growth  of  the  second-harmonic  signal  and  fit  its 
dependence  to  a  quadratic  function  of  time,  as  shown  in 
Fig.  1.  We  fit  the  data  up  to  the  point  where  it  began  to 
saturate,  i.e.  where  the  growth  was  no  longer  quadratic  in 
lime,  which  could  take  longer  than  40  minutes  at  low 
infrared  and  green  intensities.  We  performed  fits  using  two 
different  fitting  functions: 
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Fig.  1.  Quadratic  gnmih  nl'  iho  second  harmonic  signal  during 
seeding.  TTie  solid  line  is  a  til  using  Eq  The  nisei  shoves 
the  same  data  plotted  logdog  with  the  best  lit  living  Eq.  1 
(dashed)  and  using  Eq.  3  i solid i. 

=  rr  +  olTscl  or  (2) 

jpr^dued  ^  r[,  _  olTset.  (3) 

where  “offset”  is  a  small  DC  voltage  produced  by  our 
detection  system.  We  found  that  Eq.  (3)  fit  the  data  better 
than  Eq.  (2),  especially  at  early  times,  as  shown  in  the 
inset  of  Fig.  1.  We  derived  Eq.  (3)  by  assuming  an 
initially  unoccupied  level  must  first  be  populated  (with  a 
time  constant  t)  before  the  internal  DC  electric  field  can 
begin  to  form.  This  idea  was  inspired  by  recent 
suggestions  that  metasuible  levels  may  be  involved  in  the 
induced  second-harmonic  generation  effect. We  found 
experimentally  that  t  varied  from  3  minutes  at  low 
intensities  to  less  than  a  .second  at  high  intensities. 

For  all  feasible  power  levels,  we  determined  the  growth 
rate  F  and  then  fit  it  to  a  power  law  in  the  optical 
intensities:  either  Al^'  for  the  infrared  or  Bl2oj^  for  the 
green.  Figure  2  shows  the  dramatic  variation  of  F  with  the 
infrared  seeding  intensity  and  the  less  dramatic  variation  of 
F  with  the  green  seeding  intensity.  Figure  3  shows  the 
results  of  accumulating  many  such  power  law 
dependencies;  we  plot  the  c.xponent  .x  in  F  =  AI,„'  as  a 
function  of  the  green  seeding  intensity  in  Fig.  3a,  and  the 
exponent  y  in  F  =  BFio*  vs.  the  infrared  seeding  intensity 
in  Fig.  3b.  We  found  that  the  exponent  y  varied  from 
y=3.5  at  low  infrared  intensity  to  y=2.5  at  high  infrared 
intensity.  In  contrast,  the  exponent  x  varied  from  x=12.5 
at  low  green  seeding  intensity  to  x=8  at  high  green  seeding 
intensity. 

How  do  our  data  compare  with  the  current  theories  of 
second-harmonic  generation?  The  models  of  Dianov  et  al  ^ 
and  Anderson  et  al.^  postulate  a  spatially  oscillating 
photogalvanic  (jpg)  current  to  explain  the  generation  of  the 
DC  electric  field.  At  early  times  the  current  flows 
unchecked  and  the  DC  electric  field  grows  linearly  in  time: 
Eje  «  jpg  t.  In  our  experiments  we  measured  the  second- 
harmonic  signal  power,  which  is  proportional  to  Ix'^'l’of 
Eq.  (1),  so  that  according  to  the  photogalvanic  current 
models*  ’,  we  measured  Ij^^l^.  The  predicted  intensity 

dependencies  for  our  parameter  F  from  the  model  of  ref.  (8| 
is  then; 

r  «  {/2.  1 1 }  il  i.-J) 


while  from  the  model  of  ref.  [dj  n  is: 

r  -  {/:„  ll  I  ni  s-  1}  /;.  (5) 

In  Eqs.  t4)  and  (5)  the  terms  inside  the  braces  are 
proportional  to  IjpgF,  while  the  reading  process  itself 

imparts  the  additional  quadratic  variation  to  the 
measured  growth  rate  F.  This  exoa  quadratic  variation  is 
not  normali/ed  out  of  the  data  displayed  in  Figs.  2  and  3. 
If  we  include  this  built-in  quadratic  dependence  on  1(„,  then 
the  anisotropic  photoioni/.ation  model’  predicts  a  power 
law  F  oc  with  the  exponent  ranging  from  x  =  4-8. 
.Another  model*  gives  x  =  4.  The  trend  in  Fig.  3  is 
consistent  with  the  model  of  Anderson  et  al  ,  in  that  the 
power  law  exponent  drops  to  a  lower  xalue  as  the  green 
seeding  intensity  is  increased,  However,  even  at  the 
highest  obtainable  green  intensity  we  still  observe  x=H. 
This  data  clearly  indicates  a  higher-order  dependence  on 
than  predicted  by  either  model. 

There  have  been  previous  reports”’’-*  of  the  intensity 
dependence  of  seeded  second-harmonic  generation  in 
germanium-doped  fused  silica.  Kamal  et  al.'‘^  performed 
experiments  quite  similar  to  the  ones  reported  here, 
however  their  experiments  were  carried  out  in  a  specially 
grown  fiber.  If  their  data  are  analyzed  using  the  growth  rate 
parameter  F  discussed  here,  then  they  observed  x  =  4  for 
the  dependence  of  the  rate  on  infrared  intensity,  and  y  =  3 
for  the  dependence  on  green  intensity.  Their  observed 
dependence  of  the  rate  on  the  green  intensity  agrees  with 
ours,  but  their  infrared  dependence  is  quite  different. 

Lawandy  and  Selker'®  reported  a  distinct  infrared 
intensity  threshold  below  which  no  second-harmonic 
generation  was  observed  in  samples  that  had  been  seeded  for 
only  20  minutes.  Our  data  also  shows  a  very  dramatic 
dependence 


Fig.  2.  (a)  Dependence  of  ihe  rate  F  on  the  infrared  seeding 
power,  with  the  green  seeding  power  held  fixed  at  6  mW.  (b) 
The  dependence  of  ihe  rate  F  vs.  the  green  seeding  power,  with 
the  infrared  power  held  constant  at  3  W. 
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Fig.  3.  (a)  Dopendi’ncc  of  ihe  infrared  power  law  exponent  "x” 
on  the  green  seeding  set-point  intensity,  (b)  Dependence  of 
the  green  power  law  exptinent  “y”  on  the  infrared  set-point 
intensity. 

of  the  sccond-htirmonic  ciTcct  with  the  inlrtircd  inicn.sity, 
but  at  intensities  100  times  greater  than  the  threshold  cited 
in  ref.  [15|.  Rather  than  observe  a  threshold  intensity 
below  which  no  effect  occurs,  we  find  that  at  low  infrared 
intensities  the  formation  rate  of  the  DC  electric  field 
simply  becomes  extremely  slow. 

Our  observed  green  intensity  scaling  daui  are  not  in  gross 
disagreement  with  the  anisotropic  phoioioni/ation  model  of 
Anderson  ei  al^  However,  the  inlhired  scaling  daia  requires 
more  infrared  photons  than  this  model  postulates.  We 
examined  our  experiment  to  see  if  systematic  errors  caused 
the  higher-than-cxpccted  1(„  variation.  First  we  checked 
whether  the  beam  spot  si/c  entering  the  focusing  lens 
changed  as  we  adjusted  the  infrared  power.  We  measured 
the  infrared  beam  profile  and  found  no  systematic  variation 
of  the  beam  width  to  within  experimental  error  (0.6%). 
Additiortally,  we  checked  to  see  if  self-focusing  or  thermal 
Icnsing  changed  the  focal  spot  size  as  the  power  was  varied 
and  we  found  no  variation.  We  intentionally  displaced  the 
infrared  beam  from  the  green  seeding  beam  by  ~6  |im  at  the 
focus  to  determine  whether  thermal  beam  steering 
contributed  to  the  intensity  dependence,  and  found  that 
although  the  infrared  pov/er  law  exponent  “x”  was  lower 
with  the  beams  offset,  it  was  still  within  the  error  bar  of 
Fig.  3a. 

Because  the  diameters  of  our  tightly  focused  optical 
beams  change  as  they  propagate  through  the  glass  preform, 
their  intensity  varies  in  both  the  transverse  directions  and 
along  the  propagation  direction  and  therefore  so  should  the 
growth  rate.  However,  we  performed  numerical 
calculations  and  found  that  the  measured  signal  should  still 
grow  quadratically  in  time.  Because  only  the  most  intense 
part  of  the  beam  contributes  to  the  signal  at  early  times, 
the  intensity  scaling  laws  of  the  ciu'ly  time  growth  rate  arc 
not  affected.  Another  complication  is  that  our  laser  is  Q- 
switched  and  mode-locked,  so  that  there  is  a  range  of  peak 
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inicnsiiies  present  under  the  envelope,  and  our  measured 
rate  is  an  average  over  all  of  these  peak  intensities.  This 
averaging  over  intensities  could  mask  an  even  more 
dramatic  dependence  of  the  generation  rate  on  optical 
intensity. 

In  summary,  we  measured  the  buildup  rate  of  the  second- 
harmonic  signal  in  germanium-doped  fused  silica  optical 
fiber  preforms.  The  second-harmonic  signal  rises 
quadratically  in  time.  The  rate  coefficient  increases  with 
the  intensity  of  the  infrared  seeding  beam  as  and 

w  ith  the  intensity  of  the  green  seeding  beams  as 
Both  the  I,j  and  power  laws  arc  faster  then  predicted  by 
the  model  of  Dianov  et  atf  The  variation  of  the  rate  with 
the  intensity  of  the  green  seeding  light  seems  to  agree  with 
the  model  of  .-\nderson  ei  al.'^  reasonably  well.  However, 
the  \ariation  of  the  rate  with  the  intensity  of  the  infrared 
seeding  light  is  too  rapid  to  be  explained  by  this  model. 

We  thank  Robert  M.  Pierce  for  his  experimental  assistance. 
Vince  Dominic  acknowledges  fellowship  support  from  the 
Joint  Services  Electronics  Program. 
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We  demonstrate  a  new  technique  for  displaying  the  electric  field  autocorrelation  function  of  a 
laser  pulse  in  real  time,  using  two-beam  coupling  in  a  photorefractive  crystal.  This  technique 
can  be  used  over  the  entire  visible  and  near-infrared  regions  of  the  spectrum,  even  with  weak 
laser  beams,  and  is  automatically  self-phase  matched. 


The  coherence  length  of  a  laser  pulse  can  be  measured 
by  interfering  the  pulse  with  its  own  delayed  replica. Here 
we  measure  the  coherence  length of  a  pulse  from  a  mode- 
locked  laser,  where  we  average  over  many  pulses.  As  sug¬ 
gested  by  Johnson  eial.,"  we  make  each  laser  pulse  interfere 
with  itself  in  a  photorefractive  crystal.  However,  we  improve 
upon  their  technique  by  encoding  different  time  delays 
across  a  single  laser  beam.  so  that  the  entire  field  autocorre¬ 
lation  function  can  be  determined  in  essentially  real  time  (a 
few  milliseconds ) .  Our  technique  works  over  a  wide  range  of 
wavelengths  (limited  only  by  the  spectral  response  of  the 
photorefractive  crystal)  and  with  very  weak  laser  pulses 
( —  microwatt  average  power). 

The  technique  is  outlined  in  Fig.  1.  The  incident  laser 
beam  is  split  into  two,  with  one  beam  (the  reference  beam) 
directed  into  a  photorefractive  crystal.  The  other  beam  (the 
probe  beam )  illuminates  a  diffraction  grating  at  grazing  inci¬ 
dence.  The  first-order  diffracted  beam  has  different  time  de¬ 
lays  encoded  onto  its  different  pHjrtions.  Inspection  of  Fig. 
1  shows  that  the  left  side  of  the  diffracted  probe  beam  travels 
a  shorter  distance  to  the  photorefractive  crystal  than  the 
right  side  of  the  diffracted  probe  beam,  and  so  arrives  at  the 
crystal  first.  This  diffracted  probe  beam  is  focused  by  a  lens 
into  the  photorefractive  crystal,  where  it  intersects  the  refer¬ 
ence  beam  and  forms  an  intensity  interference  pattern,  but 
only  in  those  regions  of  the  crystal  where  the  two  beams  are 
coherent  with  each  other. 

The  intensity  pattern  creates  an  ihdex-of-refraction  ho¬ 
logram  in  the  photorefractive  crystal.  The  strength  of  the 
hologram  depends  on  the  visibility  of  the  intensity  interfer¬ 
ence  fringes,  which  depends  on  the  relative  coherence  of  the 
two  interfering  beams.  Note  that  if  the  incident  laser  pulses 
are  weak,  it  may  take  many  pulses  to  build  up  the  photore¬ 
fractive  hologram.  This  hologram  couples  the  two  beams,  so 
that  the  probe  beam  experiences  gain  (or  loss,  depending  on 
the  orientation  of  the  crystal).  A  one-dimensional  array  of 
detectors  measures  the  profile  of  the  amplified  probe  beam 
and  stores  it  in  a  computer.  The  probe  beam  gain  will  vary 
with  distance  x  across  the  face  of  the  beam,  and  will  depend  on 
the  local  degree  of coherence  of  the  probe  beam  with  the  refer¬ 
ence  beam.  As  we  show  below,  the  gain  profile  of  the  ampli¬ 
fied  probe  beam  is  simply  related  to  the  electric  field  autocor¬ 
relation  function  of  the  incident  laser  beam. 

Our  coherence  length  measurement  technique  can  be 
applied  to  a  train  of  (almost)  arbitrarily  weak  laser  pulses. 
The  photorefractive  hologram  integrates  the  intensity  inter¬ 
ference  pattern  over  time,  and  reaches  a  steady-state  ampli¬ 


tude  that  is  independent  of  the  total  intensity  of  the  incident 
beams.  ( However,  the  time  required  to  reach  the  steady  state 
will  increase  if  the  laser  intensity  is  decreased. )  For  example, 
in  the  experiments  reported  here,  the  time-average  intensity 
of  the  reference  beam  was  3.3  W/cm'  ( power,  74  m  W;  beam 
diameter,  1.7  mm).  At  this  optical  intensity,  the  response 
time  of  the  crystal  was  about  one-tenth  of  a  second,  enabling 
an  essentially  “real  time"  measurement  of  the  coherence 
length.  However,  we  also  tried  decreasing  the  incident  laser 
intensity  by  a  factor  of  1 000,  so  that  the  photorefractive  crys¬ 
tal  needed  a  minute  or  so  to  build  up  its  hologram,  and  we 
still  obtained  the  same  value  for  the  average  coherence 
length  of  each  laser  pulse. 

The  device  is  calibrated  as  follows.  The  reference  beam 
is  delayed  relative  to  the  probe  beam  by  an  optical  delay  line 
consisting  of  a  retroreflecting  prism  mounted  on  a  calibrated 
translation  stage.  We  repeatedly  take  data  with  different 
known  time  delays,  and  we  measure  the  translation  of  the 
output  beam’s  profile  across  the  reticon  array,  as  shown  in 
Fig.  2.  This  plot  is  linear,  as  expected,  and  its  slope  yields  the 
calibration  factor  of  our  apparatus. 

Here  are  some  experimental  details:  The  photorefrac¬ 
tive  crystal  is  a  single-domain  crystal  of  BaTiOs ,  which  we 
use  at  room  temperature  without  an  applied  electric  field. 
Both  the  reference  beam  and  the  probe  beam  are  polarized 


Delay  Line  for 
Calibration 


FIG.  I.  Optical  setup  for  measuring  the  average  coherence  leagth  of  a  senes 
of  optical  pulses.  The  diffraction  grating  encodes  different  tone  delays  onto 
different  pontons  of  the  probe  beam. 
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Position  on  Rclicon  (moi) 

FIG.  2.  Calibrating  the  apparatus.  The  bottom  ponkm  of  the  figure  shows 
the  output  beam  profile  translation  across  the  reticoa  as  the  time  delay  is 
changed-  The  top  graph  shows  the  linear  relationship  between  the  known 
time  delay  and  the  displacement  of  the  beam  profile. 


perpendicular  to  the  plane  of  incidence,  and  both  beams  are 
ordinary  rays  in  the  BaTiO,  crystal.  (Although  this  choice 
of  polarization  has  the  drawback  of  reducing  the  overall  two- 
beam  coupling  gain  of  the  probe  beam,  it  has  the  advantage 
of  minimizing  two-beam  se//-couplmg  between  the  various 
angular  components  of  the  probe  beam.)  An/=  25  cm  lens 
focuses  light  from  the  diffraction  grating  into  the  BaTiOj 
crystal,  which  is  placed  one  focal  length  from  the  lens.  The 
diffraction  grating  is  5.8  cm  square  with  1800  lines/mm.  The 
focused  probe  beam  has  a  full  cone  angle  of  12*  outside  the 
crystal.  The  unfocused  reference  beam  (intensity  full  width 
at  half  maximum  =  1.7  mm)  and  the  probe  beam  are  sym¬ 
metrically  incident  on  the  crystafat  exterior  angles  of  +  19*. 
The  BaTiO  1  crystal  is  oriented  with  its  c  axis  parallel  to  the 
wave  vector  of  the  photorefractive  hologram,  with  the  crys¬ 
tal’s  positive  c-axis  direction  chosen  to  amplify  the  probe 
beam.  With  this  geometry  the  maximum  intensity  gain  of  the 
probe  beam  is  2-5,  depending  on  which  BaTiO,  crystal  we 
use.  Depletion  of  the  reference  beam  was  minimized  by  using 
a  weak  probe  beam  ( power  =  5  /r W ) .  We  prevent  scattered 
light  from  striking  the  detector  array  by  inserting  slits  and 
apertures  in  both  the  reference  and  the  probe  beams.  We 
reduce  speckle  noise  by  inserting  a  rotating  plastic  difhiser 
between  the  crystal  and  the  detector  array,  and  then  imaging 
the  diffuser  onto  the  array.  Some  of  oar  BaTiO,  crystals 
have  dark  storage  times  of  over  an  hour,  so  in  principle  we 
can  integrate  laser  pulses  for  that  length  of  time,  if  desired. 
However,  in  practice  the  useful  integratioa  time  is  limited  to 
no  more  than  a  few  minutes,  due  to  slow  drifts  in  the  position 
of  the  mirrors  on  the  optical  table. 

If  the  diffraction  grating  were  replaced  by  a  mirror,  so 
that  each  portion  x  of  the  probe  beam  had  the  same  delay  r 
with  respect  to  the  reference  beam  then  the  intensity 


/oui  (x)  of  the  amplified  probe  after  two-wave  mixing  in  a 
crystal  of  length  L  would  be’' 

[/„„.(x)  -/,„(x)l//.,{x)  =  1],  (1) 

where  G(  r)  is  the  normalized  electric  field  autocorrelation 
function  of  the  incident  beam,  T  is  the  imensity  two-wave 
mixing  gain  per  cm  for  the  case  of  perfectly  coherent  cw  light 
beams,  and  a  is  the  intensity  absorption  coefficient.  How¬ 
ever,  unlike  a  mirror,  the  diffraction  grating  encodes  differ¬ 
ent  time  delays  r  on  each  portion  x  of  the  probe  beam  accord¬ 
ing  to 

r(x)  =  x(sin  -  sin  )/ccos5o^,  (2) 

where  x  is  the  transverse  distance  across  the  diffracted  probe 
beam,  c  is  the  speed  of  light  in  air,  is  the  angle  between  the 
incident  beam  and  the  grating  normal,  and  0^^  is  the  angle 
between  the  diffracted  beam  and  the  grating  normal.  At 
grazing  incidence  =  90*  and  6^,  =  CT,  so  that  Eq.  (2) 
simplifies  to  r(x)  =  x/c.  The  maximum  delay  ris  limited  by 
the  size  of  the  diffraction  grating,  and  was  ~  190  ps  with  our 
5.8-cm-wide  grating. 

Equation  ( 1 )  requires  a  measurement  of  the  incident 
probe  beam  profile  /,„  (x)  at  the  entrance  of  the  crystal.  In 
practice  it  is  simpler  to  measure  the  transmitted  probe  beam 
after  the  crystal  but  in  the  absence  of  any  two-beam  cou¬ 
pling,  coupling  (Jc).  by  blocking  the  reference  beam.  Equa¬ 
tion  ( 1 )  becomes 

=  iC(r)|  =  [e^^  _  _  j. 

^na  coupling  ( •T  ) 

(3) 

This  equation  includes  the  effects  of  Fitsnel  reflections 
at  the  crystal  faces,  which  happen  to  cancel  oat  nicely.  Note 
that  in  the  absence  of  linear  absorption  a,  Eq.  ( 3 )  becomes 
identical  to  Eq.  ( 1 ).  For  the  case  of  cohcrait  cw  illumina¬ 
tion,  |G  I  =  1,  and  Eq.  (3)  reverts  to  the  usual  expression 
describing  photorefractive  gain. 

We  tested  our  device  by  measuring  the  average  coher¬ 
ence  length  of  the  frequency-doubled  (i  =  532  nm )  pulses 
from  a  mode-locked  cw  Nd:YAG  laser  ((Coherent  Antares 
76-s ) .  This  laser  is  actively  mode  locked  and  emks  pulses  at  a 
rateof76  MHz.  Figure  3  shows  the  relative  gain  [defined  as 
the  left-hand  side  of  Eq.  ( 3 )  ]  across  the  probe  beam  for  a 
series  of  ~  10*  laser  pulses  (i.e.,  we  averaged  for  ~  1  s).  The 
solid  line  is  a  least-squares  fit  ofEq.  (3 )  oblaiaed  by  varying 
three  parameters;  the  maiifniim  two-beam  coopling  gain  F 
of  the  crystal,  a  vertical  offset  rebted  to  afaaorption,  and  the 
coherence  time  r,  of  the  incident  laser  pobe.  Here  we  as¬ 
sume  that  the  laser  pulses  are  Gaussian  in  lime,  and  we  take 
the  coherence  time  r,  to  be  the  half  width  at  half  maximum 
of  the  electric  field  autocorrelation  funetkm.*  From  this  fit 
we  obtain  a  value  of  r,  =  48  ps.  However,  we  found  that 
could  drift  slowly  with  time  by  —8  ps  over  JOmin,  which  we 
attribute  to  a  mode-locking  instability  of  our  hser. 

In  order  to  test  Eq.  (3),  wc  tried  usmg  two  different 
single-domain  crystals  of  BaTiO,  having  dbsimilar  two- 
beam  coupling  gain  coefficients  T.  We  obtained  the  same 
value  for  the  average  coherence  length  of  oar  laser  pulses, 
within  experimental  error,  using  either  crymaL  We  also  in¬ 
spected  our  laser  pulses  with  a  1/4-in. -thick  solid  glass 
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no.  3.  Plot  of  the  relative  intensity  gain  across  the  probe  beam.  Each  por¬ 
tion  of  the  beam  espenences  a  different  gam  because  it  has  been  delayed  by  a 
different  amount  of  time  relative  to  the  reference  beam.  The  horizontal  axis 
shows  the  relative  delay  of  each  portion  of  the  probe  be.>m.  The  solid  line  is  a 
best  fit  to  Eq.  (3).  and  yields  a  coherence  time  of  =  48  ps. 


Fabry-Perot  etalon;  the  appearance  of  the  ring  pattern  was 
consistent  with  the  coherence  length  measured  by  our  two- 
beam  coupling  technique.  However,  the  Fabry-Perot  also 
revealed  that  the  center  frequency  of  the  laser  jittered  on  the 
time  scale  of  a  few  seconds,  possibly  due  to  thermally  in¬ 
duced  fluctuations  in  the  laser's  effective  cavity  length. 

.  What  is  the  resolution  limit  of  this  device?  Dispersion  of 
the  grating  will  cause  spatial  broadening  of  the  measured 
output  beam  profile  For  the  geometry  and  coher¬ 

ence  times  used  here,  this  broadening  is  negligible,  but  it 
places  a  lower  limit  of  0. 1  ps  on  the  shortest  pulse  we  can 
measure.  Also,  the  finite  crossing  angle  of  the  two  beams  will 
cause  a  broadening  of  the  output  beam  profile.  The  worst 
case  would  be  for  a  crossing  angle  of  90*,  which  would  cause 
a  smeanng  in  time  equal  to  the  time  it  takes  light  to  traverse 
the  wider  of  the  two  beams  in  the  crystal.  In  our  experiments, 
we  estimate  this  smearing  to  be  ~  1  ps.  Our  calculations 
indicate  that  the  effect  of  a  finite  recombination  lime  of  the 
photorefractive  crystal  is  to  alter  the  overall  two-beam  cou¬ 
pling  gain  coefficient,  but  without  altering  the  shape  of  the 
output  profile  of  the  amplified  beam. 


In  summary,  we  have  demonstrated  a  simple  method  for 
measuring  the  coherence  length  of  a  pulsed  laser,  and  espe¬ 
cially  for  a  laser  that  produces  a  train  of  very  weak  optical 
pulses.  This  technique  does  not  requiie  phase  matching,  is 
simple  to  align,  and  can  be  used  over  the  entire  visible  and 
near-infrared  regions  of  the  spectrum.  The  spectral  response 
of  our  BaTiO,  crystals  extends  only  to  the  near  infrared,  and 
so  is  too  small  to  be  useful  at  1.06 /tm  Nd:YAG  laser  wave¬ 
length.  We  are  pursuing  the  use  of  other  photorefi^tive 
crystals  to  extend  this  pulse-measuring  technique  further 
into  the  infrared. 
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We  describe  how  optical  beams  from  a  mode-locked  laser  will  couple  in  a  photorefractive  cryital.  We  show  that 
the  two-beam-coupling  gain  coefficient  is  proportional  to  the  square  of  the  dsciric  field  correlation  of  the  in¬ 
coming  light  pulses.  Consequently  we  show  that  two-beam-coupling  expenmeots  can  measure  the  average  co¬ 
herence  length  of  mode-locked  laser  pulses.  We  also  describe  how  the  tempnrai  envelopes  of  the  mode-locked 
optical  pulses  distort  as  they  couple  and  propagate  through  the  photorefractive ayatal,  and  we  give  examples  of 
how  pulses  can  be  shaped  by  using  photorefractive  coupling. 


1.  INTRODUCTION 

The  strong  coupling  between  optical  beams  in  a  pho¬ 
torefractive  crystal  has  been  extensively  studied  for 
continuous-wave  optical  beams'  and  for  single,  short, 
high-irradiance  optical  pulses.^'*  Here  we  consider  the 
coupling  between  trains  of  mode-locked  pulses  in  a  photo- 
refractive  crystal.  We  are  especially  interested  in  photo¬ 
refractive  crystals  that  have  a  long  dark-storage  time, 
which  permits  the  light-induced  changes  caused  by  suc¬ 
cessive  optical  pulses  to  accumulate  in  the  crystal.  ^ 
consider  the  case  in  which  each  optical  pulse  is  quite  weak, 
so  that  it  takes  a  large  number  of  pulses  ( >  10*)  to  build  up 
a  quasi-steady-state  refractive-index  grating  in  the  crys¬ 
tal.  We  derive  the  time-averaged  response  of  the  crystal 
to  these  pulses.  We  show  that,  once  formed,  the  photo¬ 
refractive  grating  can  alter  the  temporal  shape  of  the 
laser  pulses  as  they  traverse  the  crystal.  We  also  show 
that  a  photorefractive  crystal  can  be  used  to  measure  the 
average  coherence  length  of  a  train  of  mode-locked  laser 
pulses,  as  was  recently  demonstrated  in  experiments  by 
Johnson  et  al.^~  and  Dominic  et  al.* 

2.  SPACE-CHARGE  ELECTRIC  FIELD 

Consider  an  infinite  and  periodic  train  of  laser  pulses, 
such  as  from  a  cw-pumped  mode-locked  laser.  Let  two 
such  pulse  trains  intersect  in  a  photorefractive  crystaL  If 
in  some  region  of  the  crystal  the  two  beams  are  coherent 
(and  so  produce  an  intensity-interference  pattern  there), 
then  charge  carriers  in  the  crystal  will  be  rearranged  by 
the  pattern  of  bright  and  dark  fringes.  The  migration  of 
these  charge  carriers  causes  a  space-charge  field  E^r.f) 
to  grow  and  eventually  reach  a  quasi-equilibrium.  Sup¬ 
pose  that  the  energy  of  each  pulse  is  small,  so  that  many 
pulses  are  required  in  order  to  reach  the  steady  state. 
Then,  once  reaches  its  quasi-equilibrium  value,  each 
subsequent  laser  pulse  causes  E^  to  make  only  small  ex¬ 
cursions  that  always  return  it  to  its  initial  value,  as  we 
show  below. 
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The  time  depesifeace  of  this  space-charge  field  is  gov- 
emned  by  a  first-erder  dinerential  equation,  as  derived  by 
Strohkendl  et  at.*  tar  the  case  of  constant  incident  inten¬ 
sity.  However,  the  same  equation  holds  even  if  the  inten¬ 
sity  is  r hanging  in  time,  as  long  as  the  peak  intensity  is 
weak  enough  to  mtmfy  the  restriction  of  Eq.  (8)  of  Ref.  91 
For  the  case  in  which  the  recombination  rate  of  free  carri¬ 
ers  is  much  shoctar  than  the  pulse  duration,  the  space- 
charge  field  evoivea  according  to 

+  aI^t)Ejr,  t)  =  ibEi(r,  t)Et'{r,  t),  (1) 

Ot 

where  u  the  magnitude  of  the  (quasi-dc)  space-charge 
field,  £|(r,  t)  and  ()  are  the  slowly  varying  complex 
optical-field  envelopes,  and  /e(t)  -  |£i(r,()|*  +  |£i(r,r)|* 
is  proportional  to  total  intensity.  [Note  that,  in  the 
absence  of  abeorptiosi,  energy  conservation  requires  the 
intensity  to  be  a  coBStant  in  space,  which  is  why  /o(t)  has 
no  spatial  dependence.]  The  quantities  a  and  b  are  com¬ 
plicated  functions  (d*  material  parameters  and  beam  ge¬ 
ometry  but  not  of  the  optical  intensities.* 

Equation  (1)  indicates  that  the  space-charge  Held  is 
driven  by  the  interference  term  £i(r,  r)£]*(r,  0  and  is 
erased  by  the  uniform  intensity  Iq.  If  mai^  pulses  are 
required  in  order  to  reach  steady  state,  then  no  one  pulse 
can  change  EL  by  a  large  amount,  and  £«  will  re^  a 
quasi-steadf'8tata  value  that  is  determined  fay  the  balance 
of  the  driving  and  erasing  terms.  Equation  (1)  holds 
when  the  free-canrier  recombination  time  is  much  less 
than  the  duration  of  a  single  mode-locked  laser  pulse,  so 
that  the  free^arrisr  density  created  fay  light  in  the  pho¬ 
torefractive  material  closdy  follows  any  variation  of  the 
light  intensity  with  time.  If  this  is  not  the  case,  the  situ¬ 
ation  becomes  naose  complicated,  for  we  must  then  couple 
Eq.  (1)  to  another  difTeraitial  equation  that  describes  ^e 
time  evolution  of  the  freo-carrier  density.  We  discuss 
this  case  at  the  seal  of  this  section. 

Because  in  general  the  incoming  laser  pulses  are  not 
identical,  we  must  average  over  the  ensemble  of  possible 
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pulses.  The  laser  pulses  may  fluctuate  both  in  amplitude 
and  in  phase.  We  choose  to  model  the  pulses  as  having 
identical  amplitudes  but  randomly  varying  phases.  To 
our  knowledge  there  is  no  satisfactory  statistical  descrip¬ 
tion  of  phase  fluctuations  of  mode-locked  laser  pulses,  so 
we  will  not  specify  the  statistics  of  the  phase  fluctuations; 
instead,  we  simply  characterize  the  optical  electric  field  of 
the  pulses  by  a  parameter  r,,  which  specifies  the  charac¬ 
teristic  length  of  time  (averaged  over  many  pulses)  during 
which  the  optical  electric  fields  evolve  sinusoidally  in  time 
without  interruption.  (An  analogy  may  be  an  electron 
drifting  freely  through  a  solid,  with  some  average  time  be¬ 
tween  collisions.) 

Using  our  amplitude-stabilized  pulses,  we  take  the 
ensemble  average  of  Eq.  (1); 

dE^tr  t)  -  - 

- - — -  +  alo(t)E^(r,  t)  =  ibEi{r,  i)£2*(r,  t) ,  (2) 

at 

where  the  overbar  denotes  an  ensemble  average. For  the 
case  of  cw  illumination,  Eq.  (2)  can  be  solved  by  a  steady- 
state  space-charge  field  that  is  proportional  to  the 
modulation  of  the  intensity  interference  pattern.  “  For 
the  case  oi  pulsed  illumination,  we  expect  the  equilibrium 
value  of  the  space-charge  field  to  depend  in  a  similar  man¬ 
ner  also  on  the  modulation  of  the  optical  intensity  and 
therefore  on  the  relative  coherence  of  the  two  optical 
pulses  at  the  crystal.  We  postulate  that  the  steady-state 
solution  of  Eq.  (2)  is 

T’  <£i(r,t)£a*(r,f)) 

a  (|£,(r,f)|^  +  |E,(r,f)|^’ 

where  the  the  notation  { )  denotes  an  average  over  time  71 
{F{t))  =  ~  [  f(T)dr.  (4a) 

Note  that,  if  the  interfering  beams  are  cw,  then  the  time 
averages  in  Eq.  (3)  may  be  dro{^>ed,  and  we  obtain  the 
usual  dependence  of  the  steady-state  space-charge  field  on 
the  modulation  of  the  intensity- interface  pattern.  “  For 
periodic  optical  pulses,  the  postulated  dependence  of 
Eq.  (3)  may  be  viewed  as  the  dc  term  in  a  temporal  Fourier 
expansion  of  £«.  We  show  below  that  this  dc  term  is  by 
far  the  largest  term  (1^  a  factor  of  10*)  in  the  expansion, 
again  because  many  pulses  are  required  in  order  to  reach 
the  quasi-equilibrium,  so  that  no  one  pulse  can  push  £« 
far  from  its  equilibrium  value. 

Using  the  notation  of  Trebino  et  ol.“  we  can  express  the 
electric-field  envelope  of  the  optical  pulses  as 

£i(r,  t)  =  Ai(r,  ()ui(r,  t) ,  (4b) 

Ei(r,t)  =  A*(r,f  -  Ti)u^t,t  -  u),  (4c) 

where  in  Eqs.  (4)  the  functions  Ai(r,  t)  and  Ajir,  t)  are  the 
deterministic  amplitude  of  the  pulses  of  beams  1  and  2, 
respectively.  The  fimctions  u  i  (r,  t)  and  Ui(r,  t)  are  statis¬ 
tic^  factors  that  contain  information  about  the  phase 
fluctuations  of  the  two  beams;  both  of  these  functions  are 
normalized  to  unity:  juil  =  |u«|  »  1.  For  convenience 
we  have  defined  beam  2  with  a  built-in  delay  time  u,  be- 
catise  in  practice  the  two  interfering  optical  beams  are 
usually  derived  fay  splitting  one  single  beam  into  two  and 
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then  delaying  one  beam  with  respect  to  the  other  by  a 
time  Ttf .  Using  this  notation,  we  may  rewrite  the  numera¬ 
tor  of  Eq.  (3)  as 

G(Td)  =  (£i(r,f)E2’(r,ri) 

=  (Ai(r,t)A2*(r.t  -  T,^)a,(r, /)u2*(r, /  -  u))-  (5) 

When  the  phase  fluctuations  are  governed  by  jointly 
wide-sense  stationary  statistics, so  that  the  quantity 
ui(r,/)u2*(r,  t  -  Trf)  is  independent  of  time,"  we  express 
Eq.  (5)  as 

G(Ti)  =  (A,(r,t)A2*(r,t  -  (6) 

where  is  the  second-order  coherence  function, 

vriiich  depends  on  the  relative  delay  between  the  two 
pulses  (and  implicitly  on  the  parameter  r.,  the  average 
phase-coherence  time  of  a  pulse): 

y‘^’(rd)  =  Ui(r,r)a2*(r,t  -  Ti).  (7) 

Substituting  Eq.  (7)  into  Eq.  (3)  yields 

—  if>  (Ai(r,t)A2*(r./  -  T,i))y'Vri) 

t-K  - ; - (o) 

a  lo 

In  order  to  verify  our  guess  that  Eq.  (8)  is  indeed  the 
quasi-steady-state  solution  for  the  space-charge  field,  we 
numerically  integrated  Eq.  (1)  for  the  case  of  a  periodic 
train  of  pulses,  with  the  amplitude  of  each  pulse  taken  to 
be  Gaussian  in  time.  At  first  we  assumed  that  each  pulse 
was  transform  limited.  Using  a  fourth-order  Runge-^tta 
method,  we  computed  how  changed  during  each  light 
pulse.  Figures  1  and  2  show  the  deviation  of  Ek  from  the 
steady-state  value  predicted  by  Eq.  (8).  In  Fig.  1(a)  we  set 
the  relative  delay  between  the  two  writing  beams  (one 
strong  and  the  other  weak)  to  be  =  -30  psec  (so  that  the 
strong  beam  arrives  before  the  weak  beam).  In  Fig.  1(b) 
the  space-charge  field  is  seen  first  to  decay  and  then  to 
grow  during  each  laser  pulse  but  always  to  return  to  its 
initial  value.  Figure  2  shows  the  result  for  =  +70  psec 
(where  now  the  strong  beam  arrives  after  the  weak  beam). 
Here  each  laser  pulse  causes  the  space-charge  field  first  to 
decay  (slightly)  when  the  weak  beam  arrives,  then  to  grow, 
and  then  to  decay  back  to  its  steady-state  value.  Note 
that  the  vertical  scale  ia  Figs.  1(b)  and  2(b)  is  in  parts  in 
10*;  the  maximum  excursion  of  the  space-charge  Held 
horn  its  steady-state  value  (which  it  maintains  between 
laser  pulses)  is  quite  small.  When  we  picked  a  value  for 
that  was  too  small,  we  found  that  after  the  next  li^ 
pulse  £«  would  be  left  at  a  slightly  larger  value  than  when 
it  started,  until  after  many  light  pulses  it  attained  the 
steady-state  value  specified  by  Eq.  (8).  Similarly,  ^en 
we  initially  picked  £k  too  large,  then  the  next  light  pulse 
would  leave  it  glightly  smaller.  Only  with  the  value  of  £« 
given  by  Eq.  (8)  would  each  light  pulse  cause  the  space- 
cfaaige  held  to  return  to  exactly  the  same  value  that  it  had 
before  the  light  pulse. 

To  generate  Figs.  1  and  2,  we  used  material  parameters 
that  are  typical  for  photorefractive  BaTiOs  and  optical 
hfamn  with  average  beam  intensities  of  tens  of  milliwatts 
per  square  millimeter  and  pulse  durations  of  70  psec  (in¬ 
tensity  FWHM).  For  this  calculation  we  assumed  that 
(a)  the  dark  decay  rate  is  much  less  than  the  pulse  repeti- 
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space-char^  field  momentarily  deviates  fitiin  its  steady-state 
value:  It  initially  decays  and  then  recovers. 

tion  rate,  (b)  the  average  modulation  of  the  photorefrac- 
tive  grating  is  much  less  than  unity,  (c)  there  are  no 
shallow  traps,  (d)  there  is  no  spatially  uniform  dc  electric 
field  in  the  crystal,  (e)  the  spatially  periodic  free-carrier 
number  density  is  much  less  than  the  spatially  periodic 
deep-trap  number  density,^  (f)  and  the  recombination 
time  of  charges  in  the  band  is  much  leas  than  the  duraticn 
of  each  laser  pulse. 

We  also  considered  the  case  in  which  the  optical  pulses 
have  phases  that  vary  randomly  in  time,  a^  we  found 
that  the  quasi-steady-state  space-charge  field  is  still  given 
by  the  correlation  of  the  two  optical  fields,  as  predicted  fay 
Eq.  (8).  Encouraged  by  these  results,  we  then  tried  relax¬ 
ing  condition  (f)  above  by  making  the  recombination  time 
bnger  than  the  pulse  duration  by  factors  ranging  from  1 
to  140  (while  keeping  constant  the  product  of  mobility  and 
recombination  time).  We  found  that  the  space-charge 
field  still  attained  the  value  given  by  Eq.  (8)  to  within 
1  part  in  10^ 

3.  COUPLED-mVE  EQUATIONS 

As  discussed  in  Section  2,  the  interference  of  two  weak, 
infinite  trains  of  optical  pulses  eventually  produces  a 
quasi-steady-state  space-charge  field  in  a  photorefracthe 
aystal.  This  electric  field  alters  the  refractive  index  of 
the  crystal;  it  creates  a  refractive-index  grating.  Here  we 


discuss  how  the  optical  beams  are  coupled  by  this  grating. 
In  particular,  as  the  beams  traverse  the  crystal  their 
temporal  shapes  are  modified  by  their  f«wTli"g  with  the 
grating,  which  then  alters  the  grating,  wii^  then  alters 
the  beams,  which  makes  the  equations  difficult  to  solve. 
Coupled-wave  treatments  ‘  of  this  problem  usually  ignore 
any  variation  in  the  temporal  envelopes  of  the  pulses  as 
they  propagate  through  the  crystal.  Here  we  include 
these  variations. 

To  derive  the  coupled-wave  equations,  we  consider  a  ma¬ 
terial  with  no  free  currents.  Maxwefffa  equations  then 
yield 

V*«(r.t)-uo^^^^  =  0,  (9) 

where  %(r,  t)  is  the  total  real  optical  electric  field  in  the 
material  and  9(r,  t)  is  the  corresponding  real  displace¬ 
ment  vector.  In  Eq.  (9)  we  ignore  the  term  V(V  • 

Now  let  the  total  optical  field  comprise  two  light  beams 
with  the  same  nominal  frequency  {t>o,  respective  complex 
electric  fields  £i  and  £2,  wave  vectors  k|  and  ki,  and  ei- 
genpolarizations  ii  and  e2  >n  the  ciystaL  The  total  optical 
electric  field  and  displacement  vectors  aie 

«(r,r)  =  Re{[ci£i(r,/)exp(ik,  •  r) 

+  eiEiix,  f)exp(ik2  •  r)]exp(-iwo^)},  (10) 


Fig.  2.  (a)  A  strong  optical  beam  follow!  a  weak  beam  by  70  psec. 
(b)  While  these  pulses  are  present,  the  spaea-ekarge  field  grows 
and  then  decays,  but  it  never  deviates  by  evaw  1  part  in  10*  from 
its  steady-state  value. 
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and  with  time: 


Yao  et  al . 


9»(r./)  =  ReP)i(r,t)exp(jki  -  r) 

+  D2(r,/)exp(ik,  •  r)]exp(-i<Dot)},  (11) 
D,(r,t)  s  eoKr)  •  i,Ei{r,t),  i  =  1,2.  (12) 


In  Eqs.  (10)-(12),  co  is  the  permittivity  of  free  space,  e,(r) 
is  the  relative  dielectric  tenaor  of  the  medium,  and  D.(r,  t) 
and  e,E,{r,t)  are  the  slowly  varying  envelopes  of  the  dis¬ 
placement  and  the  electric-fiekl  vectors  at  the  optical  fre¬ 
quency  <uo-  Note  that  Eq.  (12)  ignores  any  material 
dispersion.  For  laser  pulses  that  are  picoseconds  or 
longer,  this  approximation  is  justified  by  the  short  length 
of  typical  photorefractive  crystals.  For  example,  a  2-mm 
sample  of  barium  titanate  would  require  a  pulse  duration 
of  less  than  100  fsec  before  group-velocity  dispersion 
would  become  noticeable. 

The  interference  pattern  of  the  two  light  beams  will 
have  a  periodicity  C08(k,  •  rX  where  k,  =  ki  -  kj.  The 
dielectric  tensor  of  the  crystal  will  be  ^tered  by  the  light 
pattern  at  the  same  spatial  frequency  k, : 


?{r)  =  tr  +  Re(A?r  •  cxp{ik,  •  r)],  (13) 

where  e,  is  the  crystaTs  average  dielectric  tensor.  The 
spatially  periodic  part  of  the  dielectric  tensor  Re[A  e,  • 
exp(ik,  ■  r )]  will  cotqde  tihe  two  optical  waves.  Substitut¬ 
ing  Eqs.  (10)-(13)  into  Eq.  (9)  yields  equations  that  couple 
the  slowly  varying  electrk-field  envelopes  £i  and  Ei  of  the 
two  optical  beams: 


5  V  dt  / 

\  dz  V  dt 


E^Z,t)  = 


4aiC  cos  6 


4n]C  cos  8 


(eV- A?,-e2)£j(2,t), 
(14a) 

(ea*- Af/*e^)£i(2,t). 

(14b) 


In  Eqs.  (14),  z  is  the  distance  along  the  bisector  of  ki  and 
kj  in  the  crystal,  8  is  the  half-angle  between  ki  and  kz, 
and  i;  is  defined  as 


V  s  (e/n)co8  8. 


(15) 


For  cw  beams  the  time  derivatives  in  the  coupled-wave 
equations  (14)  are  zero,  but  for  pulsed  laser  beams  these 
time  derivatives  must  be  indoded.  Smirl  et  ai.‘  included 
these  time  derivatives  in  their  analysis  of  beam  coupling  in 
photorefractive  crystals  with  single-shot,  high-irradiance 
laser  pulses.  However,  they  did  not  discuss  the  evolution 
of  the  laser  pulses’  temporal  shapes,  and  with  good  reason: 
in  the  case  that  they  ronawiwed,  the  energy  coupling  be¬ 
tween  the  two  beams  was  ao  small  that  its  effect  on  the 
temporal  shape  of  the  laser  pulses  could  safely  be  ne¬ 
glected.  Here  the  energy  twiqiling  between  the  beams 
can  be  quite  strong,  and  it  can  drastically  alter  the  tempo¬ 
ral  shapes  of  the  li^t  pulses.  We  discuss  such  pulse  shap¬ 
ing  in  Section  6. 

In  deriving  the  coupled-wave  equations  (14),  we  made 
the  following  aasunqAiona:  The  optical  electric-field  en¬ 
velope  changes  slovdy  with  distance: 


dr* 


<c 


dz 


(16) 


a*(  r  •  e,E.) 

«  it) 

die  ■  e,£,) 

at* 

dt 

a(A£.  •  e,E,) 

dt 

«  oiiAe,  ■  e,E,\. 

(17a) 

(17b) 


We  assume  that  the  pertiuhation  of  the  dielectric  tensor  is 
small: 


Ae,  «  £r,  (18) 

and  that  all  the  frequency  components  of  the  laser  pulses 
are  perfectly  Bragg  matched  to  the  grating.  We  also  as¬ 
sume  that  these  frequency  components  all  write  the  same 
grating.  (For  a  millimeter-sized  grating,  this  implies  a 
minimum  laser  pulse  width  of  a  few  tens  of  picoseconds.) 


4.  PHOTOREFRACTIVE  GRATING 

The  space-charge  field  E^  of  Elq.  (8)  will  induce  a  change  in 
the  dielectric  tensor  by  the  electro-optic  (Pockels)  effect^: 

A?,  =  -C,  •  (R‘  kf)  •  trEu:,  (19) 

where  R  is  the  electro-optic  tensor  of  the  crystal. 
Equation  (19)  ignores  any  contribution  of  Ae,  from  elec¬ 
trons  in  the  conduction  band,  i.e.,  the  free-electron 
grating.*'^  Because  this  transient  grating  disappears  in 
the  time  between  each  laser  pulse,  at  the  intensities  used 
in  our  calculations  the  free-carrier  grating  will  be  negligi¬ 
ble  compared  with  the  steady-state  photorefractive  grat¬ 
ing,  which  accumulates  over  many  laser  pulses.  For 
example,  even  with  a  pulse  energy  as  large  as  6  mJ/cm*, 
the  peak  diffraction  efHciency  of  just  the  transient  free- 
carrier  grating  in  BaTiOs  is  o^y  10*^,’  which  is  200  times 
smaller  than  the  steady-state  diffraction  efficiency  with 
cw  beams  in  the  same  geometry.  Because  the  diffraction 
efficiency  from  a  free-carrier  grating  depends  quadrati- 
cally  on  the  f  luence  of  the  writing  pulses,*  with  a  cw  mode- 
locked  laser  having  a  fhience  of  microjoules  per  square 
centimeter,  the  diffraction  efficiency  from  the  free- 
carrier  grating  will  be  on  the  order  of  10"*,  which  can  cer¬ 
tainly  be  ignored.  Note  that  any  contribution  to  Hr  from 
an  absorption  grating  can  be  included  in  our  coupled- 
wave  equations  by  changing  the  overall  two-beam¬ 
coupling  gain  coefficient. 

Substituting  Eqs.  (8)  and  (19)  into  the  coupled-wave 
equations  (14)  yields 

Wi(2.t)  =  nu^7^£2(2.t).  (20a) 

\dz  V  dt  I  Jo 

( ^  ^  W*,  t)  =  -vti*  ^-y^£i(2,  t) ,  (20b) 

\dz  V  dt/  Jo 


where,  with  plane  waves  incident  upon  the  crystal,  the 
field  correlation  function  G  defined  in  Eq.  (5)  is  now  per¬ 
mitted  to  vary  only  akmg  the  z  direction.  In  Eqs.  (20a) 
and  (20b)  the  coupling  waiiitantB  rj.j  are 


Viz  ~ 


4nic  cos  0  a 


[ci*  ’  Hr'  (R-  kf)  •  Hr  -  Ca] ,  (21a) 


'boetal. 
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— ~ — :  ^  •  (R  •  *f)  *  •  «i]  •  (21b) 

4n2C  cos  Ba* 

We  combine  Eqs.  (20),  perform  an  ensemble  average  over 
many  pulses,  and  average  over  the  time  T  between  optical 
pulses: 


3<|£i(z.t)|-) 

dz 

3<|£t(z.0P) 

dz 


dz 


h 

(22a) 

h 

(22b) 

=  -  n2iW». 

io 

(22c) 

propagation,  as  we  show  below.  It  is  convenient  to  sepa¬ 
rate  the  incident,  deterministic  optical-field  envelopes 
A,(0,  t)  into  a  real  amptitode  and  a  normalized  complex 
temporal  profile  J,(/)  [so  that  the  peak  value  of  T,(t)l  is 
unity]: 


A.(0,/)  =  t  =  l,2.  (27) 

This  permits  us  to  define  the  normalized  field  correlation 
function  r(Ti): 


r(T,i) 


-  r,)) 


\~d)  ■ 


(28) 


Then  substitution  into  Eq,  (25)  yields 


In  deriving  Eqs.  (22)  we  eliminate  the  time  derivatives  by 
using  the  fact  that 


(|£i(At)|^  -  (|g.(0.<)i») 
{|£.(0.<)|') 


=  [exp(guz)  -  l]|r(Ti)|^ 


(29) 


=  0  (23) 

for  any  function  F  that  has  a  period  T. 

If  the  time-averaged  intensity  of  beam  1  is  much  less 
than  that  of  beam  2,  ((|£i|^)  «  (|£j|^»,  then  the  intensity 
of  the  strong  beam  remains  constant  over  the  interaction 
length  ((IE2P)  a  fo),  and  the  above  equations  have  the 
simple  solutions: 

G{z,Td)  =  G(0,T,<)exp(i7i2z),  (24) 

(|£i(a.f)|»)  -  (|£d0,f)|^  =  |G(0.T.)lW(gu/)  -  U,  ^25) 

h 

where  the  gain  coefficient  gu  is  defined  as 


Equation  (29)  is  the  central  result  of  our  theory;  it  shows 
that  after  the  weak  beam  has  traversed  a  distance  z  in  the 
crystal,  its  intensity  gain  depends  on  the  magnitude 
squared  of  the  normalised  electric -field  correlation  func¬ 
tion  r(Ti).  For  a  given  time  delay  tj,  the  function  r(T,i)  is 
a  measure  of  the  fringe  visibility  of  the  interfering  optical 
beams.  Note  that,  fcr  the  case  of  cw  light  beams  having 
infinite  coherence  time,  (r(T,<)|  becomes  unity,  and  the  in¬ 
tensity  of  the  weak  wave  grows  exponentially  with  dis¬ 
tance,  as  expected': 

(|£i(aO|^  =  <|JSi(0,0|*)exp(gi22) .  (30) 

5.  MEASURING  COHERENCE  OF  LIGHT 
PULSE 


gi3  3  2  Rfiliju)-  (26) 

Equation  (25)  states  that  the  average  intensity  gain  of 
the  weak  beam  is  proportional  to  the  magnitude  squared 
of  the  correlation  function  of  the  two  incident  beams. 
One  factor  of  the  correlation  function  appears  because  the 
two  beams  must  first  interfere  to  create  the  refractive- 
index  grating;  the  second  factor  appears  because  each 
beam,  after  diffraction  off  the  grating,  interferes  with  the 
transmitted  fraction  of  the  ot  .sr  beam.  Our  result  is  an 
interesting  contrast  to  that  cf  Trebino  et  who  ana¬ 
lyzed  the  interaction  of  two  beams  in  a  rapidly  responding 
absorptive  medium.’'’'^*  In  that  case  the  grating  was 
newly  formed  by  each  laser  ptilse  and  decayed  to  zero  be¬ 
tween  pulses.  Consequently  the  average  intensity  of 
their  diffracted  beam  was  proportional  to  Uie  fourth-order 
coherence  function  of  the  two  interacting  beams.  In  con¬ 
trast,  in  our  case  the  grating  in  the  photorefractive  crystal 
is  formed  by  a  large  number  of  pulses  and  reaches  a  steady 
state.  The  early  pulses  create  the  grating;  the  later 
pulses  diffract  off  the  grating.  Consequently,  the  inten¬ 
sity  of  the  transmitted  beam  is  proportional  to  the  product 
of  two  second-order  coherence  functions  and  not  to  the 
fourth-order  coherence  function. 

Equation  (25)  also  shows  that  the  average  intensity  of 
the  weak  beam  grows  at  a  rate  determined  by  the  electric- 
field  correlation  function  G(0,Ti)  evaluated  at  the  en¬ 
trance  face  of  the  crystal,  even  though  the  temporal 
profiles  of  the  two  beams  can  change  dramatically  during 


Now  we  show  how  a  slow  photorefractive  crystal  can  be 
used  to  measure  the  autocorrelation  of  a  fast  laser  pulse, 
as  experimentally  demonstrated  by  Johnson  et  al.^'  and 
by  Dominic  et  al.* 

Let  the  incident  laser  beam  be  split  into  two  beams  and 
recombined  in  a  photorefractive  crystal,  with  a  relative 
time  delay  between  the  two  beams.  Measure  the  cou¬ 
pling  gain  [defined  ly  Eq.  (29)]  as  a  functon  of  7^.  Ac¬ 
cording  to  Eq.  (29),  the  multing  plot  will  be  proportional 
to  the  square  of  the  electnc-field  autocorrelation  of  the 
incoming  laser  pulses. 

As  an  example,  consider  transform-limited  laser  pulses 
whose  intensity  profiles  are  Gaussian  in  time  with  a 
FWHM  of  Tf.  A  jdot  of  the  measured  intensity-coupling 
gain  versus  for  these  pulses  will  be  a  Gaussian  curve 
with  a  FWHM  of2v 

One  can  also  consider  the  case  of  a  cw  beam  with  a  lim¬ 
ited  coherence  length,  for  example,  a  beam  with  a  Gaus¬ 
sian  power  q>ectrum  (and  therefore  a  Gaussian  coherence 
function,  from  the  Wiener-4Qunchin  theorem): 

y'Vrf)  =  exp(-(V/’’.*)ln  2],  (31) 

where  r,  is  the  average  phase  coherence  time.  In  that 
case  the  plot  of  beam-coupling  gain  versus  will  trace  out 
a  Gaussian  with  a  FWHM  of  2r,. 

Finally,  we  consider  the  more  general  case  in  which  the 
light  is  pulsed  and  not  transform  limited.  The  electric- 
field  envelope  of  the  optical  wave  is 
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E(0,t)  =  £oexp[Kp(0)exp[-2(t’/r/)ln  2],  (32) 


where  the  phase  9(t)  of  the  wave  may  be  varying  randomly 
with  time.  For  this  pulse  Eq.  (28)  gives  a  coherence  func¬ 
tion  r(Ta): 


r(Ta)  =  exp 


(33) 


where,  in  analogy  with  the  cw  case,  we  assume  that  the 
phase  correL^ion  is  also  Gaussian,  so  that  the  ensemble 
average  of  the  correlation  of  the  phase  functions  [as  de¬ 
scribed  in  Eq.  (7)]  is 

y'“'(Ta)  =  exp(-(TaVr.*)ln  2],  ’  (34) 


vhere  r,  is  the  average  phase  coherence  time.  For  this 
case  of  non-transform-limited  pulses  the  measured  gain 
will  be  a  Gaussian  function  of  the  rebtive  time  delay  Ti 
with  a  FWHM“  of 


FWHM  = 


2t,t, 


(r/  +  T.^)^ 


(35) 


6.  PULSE  SHAPING 

The  photorefractive  grating  not  only  alters  the  intensity 
of  each  beam  but  also  alters  the  temporal  shape  of  the 
pulses  as  they  traverse  the  crystal  Tliese  pulse-shape 
changes  are  not  limited  by  the  response  time  of  the  crystal. 

Consider  the  case  in  which  two  optical  beams  have  al¬ 
ready  built  up  a  steady-state  grating  in  the  crystal.  Each 
beam  diffracts  off  this  grating  and  interferes  with  the 
other  beam,  which  drastically  alters  tbe  temporal  pulse 
shapes  of  both  beams,  as  shown  in  Pig.l 

To  describe  pulse  shaping  by  photore&nctive  two-beam 
coupling  quantitatively,  we  use  a  Fourier  rq>re8entation  of 
the  optical-pulse  envelopes  Ex  and  Ej: 


lapet  Oatpat 

Fig.  3.  Schematic  ilhiatrating  pulse  shephig  by  taro-beam  cou¬ 
pling  in  a  pbotcrefractive  cryataL  The  equan  tamporal  profiles 
of  tbs  two  incident  pulses  are  altered  aftw  coapling  in  the  crya¬ 
taL  [A  steady-state  grating  has  already  been  bi^<q)  in  the  crys- 
taL  and  each  new  li^  pulse  barely  pertuibs  diis  grating  (see 
Figs.  I  and  2)].  Note  that  when  the  two  paisas  overlap  in  tims 
the  tranamitt^  field  of  pulse  2  eonatruetiveiy  iatarferes  with  the 
difliractsd  field  of  pulse  1,  while  the  tranamktsd  field  of  pulae  1 
destnietivehr  interferes  with  the  diffracted  field  of  pulse  2. 
Whan  the  two  pulses  do  not  overiap,  there  is  as  interference 


Ei.2{zJ)  ~  2(z,w)e'““'dtv. 

(36) 

In  the  frequency  domain  the  coupled-wave  equations  (20) 
become 

dF,(z,u>)  G{z,Td)  ^  ^ 

-  +t7i2  f'2(z,  a>), 

dz  lo 

(37a) 

dFiiz.to)  .G*(2.’’<i)  „ , 

,  -  V21  ,  Fi(z,w), 

dz  lo 

(37b) 

where  for  convenience  we  have  defined 

/  iwz  \ 

Fi(z,w)  =  expl — —  jEiiz.w), 

(38a) 

1  icoz  \ 

F2{z,tD)  =  expl j£2(2,uj). 

(38b) 

Equations  (37a)  and  (37b)  are  solved  using  Eq.  (24)  and 
the  boundary  conditions 


ai)  =  , 

(39a) 

f2(0*  t*>)  “  ^2**  * 

(39b) 

dFiiz,  u)) 

=  +7^12 

G(0,T,) 

(39c) 

dz 

T  ^3  ’ 

1*0 

h 

dFiiz,  w) 

-  ’“’721 

,G*(0,r,) 

T 

(39d) 

dz 

«*0 

/o 

and  then  transforming  the  solutions  F\(z,w)  and  £2(2,  w) 
back  into  the  time  domain.  The  resulting  solution  for  the 
eleetric-neld  envelope  of  each  pulse  as  it  propagates  a  dis¬ 
tance  z  through  the  crystal  is  then 


Ei(z,t)  =  £i|o,f  -  jco8|^(r,Trf)| 


£2(^9  0  -  EiX  0, 


>•'  -  f ) 


cos|^(z,Tj)| 


r*(r.,) 


J.(o,f-^)si 


|r(T,)| 

the  function  P(z,Td)  is  defined  as 

«*.Td)  «  ^^(e-  -  1) 


(|gi(0.0l^ 

(|£|(0.f)|’> 


r(T„)(c*^  - 1). 


(40a) 


8iniP(z,Tj)|,  (40b) 


(41) 


In  Eqp.  (40)  we  have  assumed  that  the  time-averaged  in- 
tsamty  of  beam  1  is  much  leas  than  that  of  beam  2;  i.e., 
dJfij*)  «  (|£||*).  For  simplicity  we  have  also  assumed 
tbit  both  of  ^  coupling  constants  nu  and  tju  are  real 
and  equal,  so  that  v  ^  Vit  ^  Vn*  is  the  normalized 
fifrid  correlation  of  beanu  1  and  2,  as  defined  in  Eq.  (28). 
Beam  1  will  either  gain  or  lose  energy  during  its  propaga- 
tiOB  through  the  crystal  according  to  whether  the  sign  of  17 
is  positive  or  negative. 


'&oetai. 
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Equations  (40)  show  that  Ei(z,t)  is  a  weighted  sum  of 
the  input  fields  £i[0,  f  -  (z/v)]  and  Ej[(0,  t  -  (z/v)]  at  the 
entrance  face  of  the  crystal,  with  a  similar  ezpressioo  for 
beam  2.  The  relative  weights  in  the  sum  depend  cm  the 
correlation  between  the  two  input  fields,  their  relative 
time  delay,  and  the  ratio  of  their  average  intensity.  To 
confirm  Eqs.  (40)  we  uaed  Eq.  (40a)  to  calculate  the  time- 
averaged  intensity  of  the  weak  beam  at  an  arbitrary  posi¬ 
tion  z,  and  we  obtained  the  same  result  as  with  Eq.  (25). 

As  an  example  of  how  two-beam  coupling  can  dramati¬ 
cally  alter  the  shape  of  an  optical  pulse,  consider  a  train  of 
short  pulses  (£%)  and  a  train  of  somewhat  longer  pulses 
(El)  incident  upon  the  crystal.  The  short  and  the  long 
pulses  are  assumed  to  have  the  same  spectral  width.  As- 
siune  that  the  short  pulse  is  transform  limited  and  (jsus- 
sian  in  time  and  that  the  long  pulse  is  also  Gausaian  in 
time  but  is  chirped.  Let  these  two  pulse  trains  interfere 
in  a  photorefractive  crystal.  Orient  the  c  axis  of  the  crys¬ 
tal  so  that  the  longer  pulse  loses  energy  as  it  propagates 
through  the  crystal  (i.e.,  make  the  coupling  constant 


Time  (picoseconds) 


Time  (pkoeecoods) 

Fig.  4.  PuIm  shaping  by  using  Gsuiaaian  pulses  with  the  saaie 
spectral  width.  Tlwpeaka  of  beams  rwh  the  sntraneofiKs 
of  the  crystal  at  the  same  time  (r<  >  0),  but  the  duraticB  of 
beam  1  is  10  times  grsater  than  that  of  baam  2.  (a)  Tbs  tsmpo- 
ral  profile  of  beam  1  before  (dotted  curve)  and  after  (solid  carve) 
the  photorefractive  crystal.  Note  that  beam  2  couples  snaqcT  ^ 
of  beam  1.  Ws  assume  an  incident  average  intensity  raiiD  of 
ItHi  >  100  and  a  coupling  strength  pL  >  -1.23.  (b)  The  tampo- 
rsil  profile  of  beam  1  befmre  (dotted  curve)  and  aftm  (solid  esrve) 
coupling  in  the  crystal.  Here  we  reverse  the  direction  of  snwgy 
coupling  in  the  crystal,  so  that  beam  2  couples  energy  into 
beam  L  Wi  use  an  incident  average  intensity  ratio  of /t//i  >  100 
and  a  couping  strength  t)L  >  ■*•1.23. 


Time  (picoseconds) 


Time  (picoseconds) 

Fig.  5.  Pulse  shaping  in  a  photorefractive  crystal  by  using  equal- 
width,  transform-limited  Gaussian  pulses  (FWHM  inten¬ 
sity  w  70  psec).  The  curves  show  the  temporal  envelope  of 
inddent  bMm  1  before  the  cryetal  (dotted  curves)  and  after  the 
eryetal  (solid  curves),  (a)  Beam  1  arrives  50  psec  before  beam  2. 
The  coupling  strength  is  set  at  pL  >  -2  (so  that  beam  1  loses 
energy),  and  the  average  intensity  ratio  ia  set  at  U/h  ^  10. 
(b)  Bmm  1  arrives  86  peec  before  bew  2.  The  coupling  strength 
is  set  at  rjL  *  -f  1.5  (so  that  beam  1  gains  energy),  and  the  aver¬ 
age  intensity  ratio  is  set  at  Ij/Ii  >  100. 


Time  (picoseconds) 

Fy.  6.  Pulse  shaping  in  a  photorefractiva  crystal  by  using  Gaus¬ 
sian  pulses  with  the  same  temporal  width  (FWHM  intensity 
70  paec).  The  temporal  envelope  of  the  incident  beam  (1)  is 
shown  before  the  crystal  (dotted  curve)  and  after  the  cryetal  (solid 
curve).  Here  beam  1  is  transform  limited,  but  beam  2  is  chirped 
and  hu  a  spectral  width  8  times  that  of  beam  1.  The  time  delay 
between  the  two  beama  ia  zero,  and  the  coupling  strength  is 
nL  >  1.5  (beam  1  gaina  energy). 


23M  J.  Opt.  Soc.  Am.  BA'oL  7,  No.  12,  December  1990 


Yao  et  al 


t)  <  0).  For  these  conditions  Fig.  4(a)  shows  that  the  long 
pulse  will  have  a  deep  hole  carved  out  of  its  center:  a 
dark  pulse.  Such  a  pulse  may  prove  useful  for  the  experi¬ 
ments  with  dark  pulses.*^’  (^e  Appendix  A  for  further 
discussion  of  this  calculation.) 

If  we  reorient  the  crystal  so  that  the  longer  pulse  gains 
energy  during  propagation  (tj  >  0),  then  it  will  emerge 
with  a  sharp  spike,  as  shown  in  Fig.  4(b).  Beam  shaping 
will  also  occur  if  the  two  pulses  have  similar  temporal  en¬ 
velopes,  as  shown  in  Figs.  5  and  6.  In  Fig.  5  we  chose  the 
two  beams  to  be  identical,  transform-limited  Gaussians. 
In  Fig.  6  we  chose  a  transform-limited  weak  beam  and  a 
strong  beam  that  is  linearly  chirped  and  with  a  somewhat 
wider  spectral  width.  Note  the  strong  ringing  caused  by 
the  temporal  overlap  of  the  two  pulses. 

7.  CONCLUSION 

In  summary,  we  have  derived  an  analytical  expression  for 
the  two-beam-coupling  gain  of  a  photorefractive  crystal 
when  it  is  illuminated  by  two  light  beams  that  contain  in¬ 
finite  trains  of  ampiitude-stahilized  pulses.  We  demon¬ 
strated  that  the  photorefractive  space-charge  field  reaches 
a  quasi-steady-state  vahie  proportional  to  the  local  fringe 
visihility  in  the  crystal  Vk  find  that  the  time-average 
intensity  gain  of  each  beam  is  proportional  to  the  mag¬ 
nitude  squared  of  the  electric-field  correlation  of  the 
interfering  beams  at  the  entrance  face  of  the  crystal. 
Consequently,  the  oorrelatian  between  two  optical  pulse 
trains  can  be  measured  fay  performing  two-beam  coupling 
experiments  in  a  pfaotarefractive  crystal.  We  have  also 
shown  that  optical  pubes  will  be  shaped  by  their  coupling 
in  a  photorefractive  crystal,  and  we  have  derived  analytic 
expressions  describing  how  each  pube  transforms  during 
propagation  through  the  crystal. 


where  t^i  is  the  FWHM  intensity  width  of  the  long  pulse 
and  do  is  a  constant  phase.  Both  pulses  are  Gaussian 
with  the  same  spectral  width  =  (4/Tp2)ln  2,  as  can  be 
shown  by  taking  their  Fourier  transforms.  Because  the 
original  pulse  is  assumed  to  be  transform  limited,  the  co¬ 
herence  function  y'*'  is  unity.  The  corresponding  correla¬ 
tion  function,  calculated  using  Eq.  (28),  becomes 


If  the  average  intensity  of  beam  1  is  much  smaller  than 
that  of  beam  2,  e.g.,  ()£\|^)  «  ((£21  ^),  then  substitution 
into  Eq.  (40a)  yields 

£,(z,/)  =  £io{9’i[t  -  (z/v)] 

+  Vqr(T,,)(e'''  -  l)J2{f  -  {zjv)  -  j,]}.  (44) 

Substituting  Eqs.  (42)  and  (43)  into  Eq.  (44),  we  can  calcu¬ 
late  the  intensity  of  pulse  1  as  a  function  of  time  and  dis¬ 
tance  in  the  crystal.  The  results  are  shown  in  Figs.  4  and 
5.  [In  Fig.  5  both  beams  are  transform  limited,  so  that 
q  =  1  in  Eqs.  (41)-(44).] 

In  Fig.  6  we  consider  pulses  with  the  same  temporal 
width  but  different  spectral  widths  and  with  a  linear 
chirp  on  the  spectrally  wider  beam.  If  the  spectral  width 
of  pulse  2  is  m  times  that  of  pulse  1,  then  the  two  pubes 
can  be  expressed  as 

J,(t)  *  exp|-2  In  2-^j  >  (45) 

~  u)  =  e*P  -2  In  2 


APPENDIX  A;  PULSE-SHAPE 
CALCULATIONS 

We  can  use  Eqs.  (40a)  and  (40b)  to  calcubte  how  the  inten¬ 
sities  of  two  interacting  pubes  will  change  with  time  be¬ 
cause  of  beam  coupling.  In  general,  we  must  first  choose 
the  shapes  of  the  incident  pubes  before  the  crystal  and 
then  calculate  the  normalued  electric-field  correlation 
function.  In  the  exangib  shown  in  Fig.  4,  we  chose  one 
pube  to  be  transform  limitwd  and  Gaussian  in  time  and 
the  other  pube  also  to  be  Gauaaian  in  time,  with  the  same 
spectral  width,  but  chirped.  We  might  generate  such 
pubes  by  splitting  a  transform-limited  Gaussian  pulse 
into  two  pubes  and  — wding  one  of  the  pubes  through  a 
Treecy  grating  pair.*  Tbu  pube  will  be  chirped  and 
temporally  stretched.  If  the  dwtion  t,i  of  the  stretched 
pube  (1)  b  q  times  the  duration  of  the  unstretched 
pube  r,i  s  qr,,,  then  the  normalized  complex  temporal 
amplitades  of  these  pubes  can  be  expressed  as 


’  -2  In  2(1  -t-  i 

l>  - 7-5 


Vg^  -  l)t* 


Pube  1  b  transform  hmited,  and  the  coherence  function 
b  unity.  The  corresponding  correlation  function,  cal¬ 
culated  using  Eq.  (28),  b 


Substituting  Eqs.  (45)-(47)  into  Eq.  (44)  and  setting 
g  *  I  (the  two  pubes  are  given  the  same  temporal  width), 
we  obtain  Fig.  6. 


ACKNOWLEDGMENTS 

We  thank  R.  Pierce,  J.-M.  C.  Jonathan,  D.  Mahgerefteh, 
and  C.  IL  Caves  for  fruitful  discussions.  This  research  is 
supported  fay  the  Joint  Services  Electronics  Program  and 
contract  F49620-8S-C-9905  of  the  U.S.  Air  Force  Office  of 
Scientifk  Research.  V  Dominic  gratefully  acknowledges 
fellovnhip  support  from  the  National  Science  Foundation 
and  the  Joint  Services  EVectronics  Program. 


YMeta/. 

REFERENCES  AND  NOTES 

1.  L.  Solymar,  “Theory  of  volume  holographie  fonnaticn  in  pho- 
torefractive  crystals.”  in  Elect ro-Optie  and  Photorefractive 
Materials.  P.  Gunter,  ed.  (Springer-VMa(,  Berlin,  1987), 
pp.  229-245. 

2.  A.  L.  Smirl.  G.  C.  Valley.  K.  M.  Bohnert.  and  T.  K  Boggeaa. 
Jr.,  “Picosecond  photorefractive  and  free- carrier  transient 
energy  transfer  in  GaAs  at  1  IEEE  <1-  Quantum  Elec¬ 
tron.  24.  289-302  (1988). 

3.  A.  L.  Smirl.  K.  Bohnert,  G.  C.  Valley,  R.  A.  Mullen,  and 
T.  F  Boggess.  “Formation,  decay,  and  ersnure  of  photorefrac¬ 
tive  gratings  written  in  barium  titanate  by  picosecond 
pulses.”  J.  Opt.  Soc.  Am.  B  6,  606-615  (19S9). 

4.  A.  L.  Smirl.  G.  C.  Valley,  R.  A.  Mullen,  K.  Bohnert, 

C.  D.  Mire,  and  T.  F.  Boggess,  “Pkoaecoad  photorefractive 
effect  in  BaTiOs,”  Opt.  Lett.  12,  501-503  (1987). 

5.  A.  M.  Johnson,  A.  M.  Glass,  W  M.  Simpenn.  R.  B.  Bylsma, 
and  0.  H.  Olson,  “Microwatt  picosecond  pulse  autocorrelator 
using  photorefractive  GaAs:Cr,”  in  OSA  Annual  Meeting, 
Vol.  11  of  1988  Optical  Society  of  America  Technical  Digest 
Series  (Optical  Society  of  America,  Washington,  D.C.,  1988), 

p.  128, 

6.  A.  M.  Johnson,  A.  M.  Glass,  W  M.  Simpson,  and  D.  H.  Olson, 
“Infrared  picosecond  pulse  diagnostics  using  photorefractive 
beam  coupling,”  in  Conference  on  Lasers  and  Electro-Optics, 
Vol.  II  of  1989  (^tical  Society  of  Anmrica  Technical  Digest 
Series  (Optical  Society  of  America,  Washington,  D.C.,  1989), 

p.  226. 

7.  A.  M.  Johnson,  W  M.  Simpson,  A.  M.  Glass,  M.  B.  Klein, 

D.  Rytz,  and  R.  Trebino,  “Infrared  picosecond  pulse  correla¬ 
tion  measurements  using  photorehiictitre  beam  coupling  and 
harmonic  generation  in  KNbOj  and  BsTiOi,”  in  OSA  Annual 
Meeting,  Vol.  18  of  1989  Optical  Society  of  America  Techni¬ 
cal  Digest  Series  (Optical  Society  of  America,  Washington, 
D.C.,  1989),  p.  53. 

8.  V  Dominic,  X.  S.  Yao,  R.  M.  Pierce,  and  J.  Feinfaerg,  ‘Measur¬ 
ing  the  coherence  length  of  mode-lockad  laser  pulses  in  real 
time,"  Appl.  Phys.  Lett.  56,  521-523  (1990). 

9.  F  P.  Strohkendl,  J.  M.  C.  Jonathan,  and  R.  W  Hellwarth, 
‘Hole-electron  competition  in  photorefractive  gratings," 
Opt.  Lett.  11,  312-314  (1986). 

10.  Here  we  use  the  fact  that  the  enaemhls  average  of  the  de¬ 
rivative  of  a  function  equals  the  dativative  of  the  ensemble 
average  of  the  function.  See,  for  example,  A.  Papoulis,  Proba¬ 
bility,  Random  Variables,  and  Stoehastie  Process  (McGraw- 
Hill,  New  York.  1965),  Chap.  9  pp.  314-31& 

11.  J.  Feinberg,  ‘Opticad  phsM  copjugatian  in  photorefractive 
materials,”  in  Optical  Phase  Conjugation.  R.  Fisher,  ed. 
(Academic,  New  York,  1983),  pp.  417-443. 

12.  R.  Trebino.  E.  K.  Gustafson,  and  A.  B.  ^iegiiiaii,  ‘Fourth- 
order  partial-coherence  effects  in  the  formation  of  integrated- 
intensity  gratings  with  pulsed  light  sources,"  J.  Opt.  Soc. 
Am.  B  3,  1295-1304  (1986). 

13.  J.  W  Goodman,  Statistical  Optics  (WQey,  Now  York,  1985), 
Chap.  3,  pp.  63-68. 

14.  For  example,  if  both  of  the  optical  fields  are  derived  from  the 
aame  mode-locked  laser,  then  Eq.  (6)  will  bold  as  kaag  as  the 
laser  exhibits  no  long-term  changas  in  the  statistics  of  its 
phase  fluctuations  during  an  experimant.  However,  if  the 


VoL  7,  No.  12/December  1990/J.  Opt.  Sot  Am.  B  2355 

alignment  of  the  mirrors  of  the  laser  should  drift  with  time, 
for  example,  then  the  laser's  pulse  statistics  might  change, 
and  Eq.  (6)  would  no  longer  follow  from  Eq.  (5). 

15.  The  coupling  coefficients  tih  and  tiq  in  tte  coupled-wave 
equations  i21a)  and  (21b)  will  be  sightly  difliwsnt  if  the  term 
V(V  ■  g)  is  retained.  Namely,  cos  9  will  be  replaced  by 
cos  9  —  Re{(z  ■  e,)(l,  ■  e,*)]  with  i  ^  1  in  Eq  (21a)  and  i  •  2 
in  Elq.  (21b),  where  Si  and  it  are  the  unit  vecSsss  along  ki  and 
k],  respectively,  i  is  a  unit  vector  in  the  z  Erection,  and  we 
assume  that  the  optical  field  amplitudes  vary  snly  along  the  z 
direction. 

16.  A.  V.  Alekseev- Popov,  A.  V.  Knyaz’kov,  and  A.  S.  Saikin, 
‘Recording  volume  amplitude-phase  bologranu  in  a  lead- 
lanthanum  zirconate-titanate  ceramic,"  Sea  Tech.  Phys.  Lett 
9,  475-477  (1983). 

17.  K.  Walsh,  T.  J.  Hall,  and  R.  E.  Burge.  ‘Infhsnce  of  polariza¬ 
tion  state  and  absorption  gratings  on  phelaefractive  two- 
wave  mixing  in  GaAs."  Opt  Lett.  12,  1026-1028  (1987). 

18.  Y.  Lee,  “Studies  of  the  photorefractive  eflart  in  barium  ti¬ 
tanate:  higher-order  spatial  harmonics  and  two-beam  en¬ 
ergy  coupling,"  (University  of  Southern  Cahfomia.  Los 
Angeles,  Calif,  1989). 

19.  R.  M.  Pierce,  R.  S.  Cudney,  G.  D.  Bacher,  tasi  J.  Feinberg, 
“Measuring  photorefractive  trap  density  wdlsMt  the  electro¬ 
optic  effect,"  Opt.  Lett.  15,  414-416  (1990);. 

20.  H.  J.  Eichler,  U.  Klein,  and  D.  Langhans,  *Coherence  time 
measurement  of  picosecond  pulses  fay  ligkt-aduced  grating 
method."  Appl.  Phys.  21,  215-219  (1980). 

21.  Z.  Vardeny  and  J.  Tauc,  ‘Picosecond  cohsrTe  coupling  in 
the  pump  and  probe  techniquq’  Opt  (^omamn.  39,  396-400 
(1981). 

22.  R.  Baltrameyunas,  Yu.  Vaitkus,  R.  DanelynB,  M.  Pyatrauskas, 
and  A.  Piskarskas,  “Applications  of  dynamic  holography  in 
determination  of  coherence  times  at  single  yanaecond  light 
pulses,"  Sov.  J.  Quantum  Electron.  12,  1252-1254  (1982). 

23.  S.  L.  Palfrey  and  T.  F.  Heinx,  ‘Cohanak  interactions  in 
pump-probe  absorption  measurements:  the  effect  of  phase 
gratings,"  J.  Opt.  Soc.  Am.  B  2,  674-678  (19Kk 

24.  W.  L.  Nigfaan,  Jr.,  T.  Gong.  L.  Liou,  and  P.  M.  Fauchet,  “Self¬ 
diffraction:  a  new  method  for  cbaracterimSiaai  of  ultrashort 
laser  pulses,"  Opt.  Commun.  69,  339-344  (1289). 

25.  In  Ref  8  we  defined  the  coherence  time  at  the  pulse  as  one- 

half  of  the  FWHM  of  the  field  correlatiim  Aa^on.  This 
definition  leads  to  a  coherence  time  »  r,  >  (2  In  Zl/fmli/)  for 
an  optical  field  having  a  Gaussian  power  igaetral  density 
with  a  FWHM  of  Av  (Hz).  With  this  definstiM  the  coherence 
time  Tc  equals  the  pulse  duration  r,  for  a  taaaaform-limited 
pulse.  A  more  common  definition**  of  r,  hr  cw  beams  is 
T,  =  / !rir(Ti)|*  dre.  In  this  case  t,  =  [(2  la  l/(Ai/)  for 

an  optical  field  having  a  (Sauasian  spectraaa  with  a  FWHM 
of  Ak  (Hz). 

26.  D.  Krokel,  N.  J.  Halas,  G.  Giuliani,  and  H  Grischkowsky, 
“Dark-pulse  propagation  m  optical  fibers,"  ItiyB.  Rev.  Lett. 
60.  29-32  (1988). 

27.  A.  M.  Weiner,  J.  P.  Heritage,  R.  J.  Hawkins,  R.  N.  Thurston. 

E.  M.  Kirschner,  D.  E.  Leaird,  and  W  J.  Taminaan.  “Experi¬ 
mental  observation  of  the  fundamental  dark  saiitoa  in  optical 
fibers,"  Phys.  Rev.  Lett.  61,  2445-2448  (Mik 

28. ;  E.  D.  Treacy,  “Optical  pulse  comproaaian  with  diffraction 

gratings,”  IEEE  J.  Quantum  Electron.  <9^454-458  (1969). 


